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The principle of corresponding states in modern form has been applied to the following 


properties: the critical state, the virial coefficient, the Boyle point, the densities of coexistent 
phases, the vapor pressure of the liquid, the entropies of evaporation and of fusion, the coeffi- 
cient of thermal expansion of the liquid, the triple-point temperature and pressure, the heat 
capacity of the liquid, and the surface tension of the liquid. It has been shown that argon, 
krypton, xenon, and with less accuracy neon, follow the principle with respect to all these 
properties. It has further been shown that nitrogen, oxygen, carbon monoxide, and methane 
follow the principle with fair accuracy as vapors and as liquids, but not as solids. The relations 
between surface tension, temperature, and densities have been analyzed empirically. For the 
“ideal’’ substances under consideration Katayama’s modification of Eétvos’ relation holds 
good, but McLeod’s relation does not; in the relation y~ (p:—p,)*, the exponent s is not 4 but 


much more nearly 33. 


1. INTRODUCTION 


HE principle of corresponding states may 

safely be regarded as the most useful by- 
product of van der Waals’ equation of state. 
Whereas this equation of state is nowadays 
recognized to be of little or no value, the principle 
of corresponding states correctly applied is ex- 
tremely useful and remarkably accurate. 

Pitzer' has stated a set of assumptions suffi- 
cient to lead to the principle of corresponding 
states and has shown that argon, krypton, and 
xenon have several properties in accordance with 
the principle. In the present paper, I shall draw 
attention to a number of further relations, not 
mentioned by Pitzer, all in accordance with the 
principle of corresponding states. I shall show 
that the principle is also followed accurately by 
minagen | and oxygen, and somewhat less accu- 


1 1K. Pitz Pitzer, a Chem. Phys. 7, 583 (1939). See also Boer 
hysica 5, 946 (1938). 
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rately by carbon monoxide and methane. For the 
sake of continuity and completeness, I shall at 
small cost of space repeat some of Pitzer’s 
contributions to the subject. 


2. ASSUMPTIONS 


I shall with the utmost brevity recall Pitzer’s 
list of assumptions and make comments on some 
of them. They are: 


I. ‘Classical statistical mechanics will be used.” 

II. ““The molecules are spherically symmetrical, either 
actually or by virtue of rapid and free rotation.” 

III. “Intramolecular vibrations will be assumed the same 
in the liquid and gas states.” 

IV. ‘‘The potential energy will be taken as a function only 
of the various intermolecular distances.” 

V. “The potential energy for a pair of molecules can be 
written Ad¢(R/Ro) where R is the intermolecular 
distance, A and R» are characteristic constants and 
¢ is a universal function.” 


The precise meaning of assumption | is ambigu- 
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ous. My interpretation of the correct condition is 
that 

(a) any distinction between Fermi-Dirac sta- 
tistics and Bose-Einstein statistics has a negli- 
gible effect ; 

(b) the effect of quantization of the transla- 

tional degrees of freedom is negligible. 
Both (a) and (b) are satisfied provided (mkT)}vt 
>h, where m denotes molecular mass, T absolute 
temperature, and v volume per molecule. As 
pointed out by Pitzer this condition excludes 
hydrogen and helium owing to their small mo- 
lecular masses and to some extent limits the 
applicability of the principle to neon. I can how- 
ever see no need to assume that degrees of 
freedom other than the translational must be 
classical. On the contrary, the vibrational degrees 
of freedom of diatomic molecules are practically 
unexcited, and yet it will be shown that such 
molecules can obey the principle of corresponding 
states. 

The effect of assumption II is that some 
diatomic and polyatomic molecules may obey the 
principle in the gaseous and liquid state, but 
cannot be expected to obey it in the solid state. 
Assumption III might be more usefully worded : 
“The intramolecular degrees of freedom are as- 
sumed to be completely independent of the 
volume per molecule.’’ Highly polar molecules 
are ruled out by conditions II, III, and IV. 
Metals and molecules capable of forming hydro- 
gen bonds are ruled out by condition IV. 

Even when conditions I to IV are fulfilled, 


there is no a priori reason why assumption V | 


should be fulfilled. In fact there is good reason to 
suppose that it does not hold accurately. The 
assumption is however accurate for large values 
of R where the intermolecular potential energy e 
is proportional to —R~*. Moreover many macro- 
scopic properties are insensitive to the precise 
form of the relation between e and R for small R. 
Consequently assumption V, though not rigor- 
ously true, turns out to be a useful approximation 
for many non-polar molecules. 


3. FORMULATION OF PRINCIPLE 


From the above assumptions it follows that the 
free energy F of a system of N molecules at the 
temperature T obeys the relation 
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where m is the molecular mass, v is the molecular 
volume, j is the partition function for the internal 
degrees of freedom of a molecule and depends 
only on the temperature 7, the principal mo- 
ments of inertia 7, and the frequencies v, of the 
normal modes. Finally ¥ is a complicated, but 
universal function of kT’/e) and v/Ro*, where Ro 
is the intermolecular distance at which the 
intermolecular energy has its minimum value 
—e. At infinite dilution y tends to unity and 
formula (3.1) becomes that of a perfect gas. 
From the thermodynamic relation 


where P is the pressure, V the volume of the 
whole system, and v the volume per molecule, it 
follows that Pv/kT is a universal function of 
kT /e9 and v/R,*. From the well known conditions 
for the critical state it follows at once that 
kT./e9 and v./Ro*, where the subscript c denotes 
critical values, are both universal constants. The 
equation of state is therefore of the form 


where u is a universal, but complicated, function. 
Insofar as the five.initial assumptions are fulfilled, 
Eq. (3.3) is valid for gas, liquid, and solid. This is 
the general expression of the principle of corre- 
sponding states. 

I shall now consider various deductions from 
formula (3.3) and examine to what extent they 
are borne out by experimental data for neon, 
argon, krypton, xenon, nitrogen, oxygen, carbon 
monoxide, and methane. 


(3.3) 


4. CRITICAL STATE 


In Table I, the first row gives the molecular 
weights, the next four rows? the critical tempera- 


? The sources of experimental data are collected in an 
appendix. . 
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Kr 


Xe Ne 


Oz co 


M 

T./°K 

p-/g 
v-/cm? mole! 
P../atmos. 
P.w./RT-. 


83.7 
209.4 


131.3 28.02 
126.0 


32.00 28.00 
154.3 133.0 
0.430 0.301 
74.5 93.2 
49.7 34.5 
0.292 0.294 


TB/°K 
Ts/T. 


~345 
2.6 


T,/°K (P =1 atmos.) 


T./°K (P=P,./50) 
T,/T-. 


81.6 

78.9 
0.613 
0.593 


L./R°K 
L./RT, 
L./RT; 


727 


1023 
9.10 
‘9.26 


10%a 


ou | 
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T,/°K 


of | om 
Co Sx 
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one 
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L;/R°K 
L;/RT, 


Dw 


Oo 
~ 


P,./atmos. 
100P./P. 


ney 


vi/em3 


v,/cm? 
v1/Us 


Cp/R at Ts 5.4 


* This is not a measured value, but an approximate value based on vapor pressure measurements. 


ture T,, the critical density p., the critical molar 
volume v,, and the critical pressure P,. The sixth 
row gives the values of P.v./RT., which van der 
Waals pointed out should have a universal value. 
It will be seen that the values for argon, krypton, 
nitrogen, oxygen, carbon monoxide, and methane 
all lie within 1.5 percent of 0.292. The value for 
neon is somewhat higher and that for xenon 
somewhat lower. In the former case the deviation 
is presumably a quantal effect ; in the latter there 
is no obvious explanation. 


5. VIRIAL COEFFICIENT AND BOYLE POINT 


The equation of state of a gas up to moderate 
pressures may be expressed in the form 


Pv/RT =1+B(T)/2, (5.1) 


where B(T) is called the second virial coefficient, 
which is positive at high and negative at low 


temperatures. The temperature Ts, at which 
B(T) changes sign is called the Boyle point. In 
the seventh row of Table I are given values of Tz 
and in the eighth values of 7,/T7.. The values for 
neon, argon, and oxygen agree within 1 percent; 
those for nitrogen, carbon monoxide, and methane 


are 5 to 8 percent lower. 


The dependence of B(T) on T is not expressible 
by any particularly simple relation. B(T)/v, is 
plotted against 7/7, in Fig. 1 for neon, argon, 
nitrogen, and oxygen, for which alone of the 
substances being considered accurate data at low 
pressures are available. It will be seen that the 
three substances are roughly fitted by the same 
curve. 


6. DENSITIES OF COEXISTENT LIQUID 
AND VAPOR 


If p: denotes the density of the liquid and p, 
that of the vapor in mutual equilibrium at the 
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Formula Ideal Ne A CHa 
20.18 39.94 P| 16.03 
1) 448 150.7 190.3 
0.484 0.5308 0.9085 1.155 0.311 0.162 
41.7 75.3 92.1 113.7 90.2 98.8 
™ 26.9 48.0 54.1 58.2 33.5 45.7 
al 0.292 0.305 0.292 0.290 0.278 0.292 0.289 
ds 121 411.5 327 491 
2.70 2.73 2.59 2.58 
he . 27.2 87.3 120.9 165.1 77.3 90.1 112.5 
10 25.2 86.9 122.0 167.9 74.1 90.1 110.5 
ut 11 0.608 0.580 0.577 0.570 0.614 0.583 0.591 
Ro 12 0.58 0.563 0.577 0.582 0.580 0.588 0.583 0.581 : 
he 13. 224* 785 1086 1520 671 820 | | 
= 14. 8.25 - 8.98 8.98 9.19 9.11 8.91 
id 15. 9.05 8.9 9.04 8.91 9.06 9.11 9.22 i= 
16. 3.85 4.91 
0.68 0.59 0.65 
18. } | 24.6 116.0 161.3 54.4 68.1 90.6 
19, 0.555 0.549 0.553 0.557 | 0.352 0.512 0.476) 
2 
) 20. 40.3 196.2 276 53.5 100.7 113 
21. 1.69 1.64 1.69 ere 2 0.98 1.48 1.25 
1e 22. 0.425 0.721 0.810 
it 23. P| 1.37 1.58 1.33 1.39 
of 34.13 42.68 
23. 29.65 37.09 
26. 1.151 1.151 
2g 27. | 5.35 5.35 6.75 6.45 7.25 
3) 
1, 
is 
n 
y 
l, 
n 
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temperature 7, while p. denotes the critical 
density, then according to the ‘principle of corre- 
sponding states one should expect p:/p. and p,/p- 
to be universal functions of 7/T7,. In Fig. 2 are 
plotted all the experimental data for the sub- 
stances under consideration. It will be seen that 
except for carbon monoxide and methane most of 
the points lie on or near a single curve. 

The curve drawn through the points for argon 
is represented by the empirical formulae 


By addition and subtraction of (6.1) and (6.2) 


one obtains 


Pit Po _ (6.3) 


2p- 

Pe 2 


Formula (6.3) expresses the well-known law of the 
rectilinear diameter. Formula (6.4) is new. The 
inaccuracy of these formulae as applied to argon 
is generally only one or two parts per thousand 
of p. or of p:. But the percentage inaccuracy in p, 
increases with decrease of temperature and be- 
comes serious below T~~0.65T7.. In fact according 
to the formulae, p, would become zero at 
T=0.5657., but this is below the triple point 
T,=0.587,. It is therefore not recommended that 


(6.4) 
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these formulae should be used for computing 
values of p,. There are however occasions when 
one requires relatively accurate values not of p, 
itself but of (p:—p,) /p-;on such occasions formula 
(6.4) in view of its extreme simplicity and 
surprisingly high accuracy has much to recom- 
mend it. An example of its use will occur in 
Section 16. 


7. VAPOR PRESSURE 


At temperatures considerably below the critical 
temperature, say T<0.657., when formula (6.2) 
for p, becomes inaccurate it is convenient to con- 
sider the equilibrium vapor pressure P rather 
than p,. According to the principle of corre- 
sponding states one should expect P/P, to be a 
universal function of 7/T.. In particular the 
temperatures 7, at which the equilibrium pres- 


sure P is one-fiftieth of the critical pressure _ 
should be corresponding temperatures for differ- 
ent substances and the ratio of T, to T. should 
have a universal value. On the other hand 7; the 
boiling points at a pressure of one atmosphere are 
not corresponding temperatures for different 
substances. In rows 9 and 10 of Table I are given 
T, the boiling point at a pressure of one atmos- 
phere, and 7, the boiling point at a pressure one- 
fiftieth the critical pressure. In rows 11 and 12 
are given the ratios 7,/T. and T,/T,. It will be 
seen that the values of the latter are, as expected, 
more nearly the same than the values of the 
former. 


8. ENTROPY OF EVAPORATION 


According to Trouton’s rule the molar entropies 
of evaporation for different substances have 


257 
100 | 
P 
95 
+ 
sot | 
4 
+ 
+ 
£0 + 
15 + Ne 
eA 
= 
v 
CO 
65h, o CH, 
al + 
Vv 
+ 
ssl Ai 
0 2 £ 10 12 14 18 22 2h 
Fic. 2. 
3) 
4) 
he 
1e 
yn. 
id 
Pog 
e- 
at 


258 


roughly the same value at their boiling-points 
under a pressure of one atmosphere. Hildebrand 
has found that the molar entropies of evaporation 
are more nearly the same at temperatures at 
which p, has the same value. Neither Trouton’s 
rule nor Hildebrand’s accords with the principle 
of corresponding states, which requires rather 
that the molar entropies of evaporation should be 
equal at corresponding temperatures. In par- 
ticular they should be equal at a standard tem- 
perature 7,=0.587, at which the vapor pres- 
sure is one-fiftieth the critical pressure. In row 13 
of Table I are given the molar heats of evapora- 
tion L, divided by the gas constant R. In rows 
14 and 15 are given the molar entropies of 
evaporation divided by the gas constant R at the 


temperatures 7, and T, respectively. It will be 


seen that the values of L./RT, are, as they should 
be, more nearly alike than the values of L./RTy. 


9. COEFFICIENT OF THERMAL EXPANSION 


An extremely sensitive test of the accuracy of 
formula (6.1) for the density p; of the liquid is 
obtained by using this formula to calculate a the 
coefficient of thermal expansion. According to 
this formula a will vary with the temperature, 
but only slowly at low temperatures. At the 
temperature 7,=0.587., the calculated value of 
a(T.—T,) is 0.287 and that of a7, is 0.68. The 
only experimental values of @ are given in row 16 
of Table I and the corresponding values of a7, 
in row 17. The agreement for nitrogen, oxygen, 
and carbon monoxide is as good as one could 
expect. In the case of argon: the agreement is 
exact. This is not an independent test of the 
principle of corresponding states, but is a demon- 
stration of the high accuracy of formula (6.1) as 
applied to argon. 


10. TRIPLE-POINT TEMPERATURE 


As pointed out in Section 2, it is only for 
monatomic molecules that the principle of corre- 
sponding states may be expected to apply to the 
solid state, because the rotational motion of 
diatomic and polyatomic molecules can hardly be 
unrestricted in a crystal. In row 18 of Table I are 
given the triple-point temperatures 7; and in 


( : —— J. Am. Chem. Soc. 37, 970 (1915) ; 40, 45 
1918). 
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row 19 the ratios of T; to T.. There is good 
agreement between the monatomic molecules of 
argon, krypton, xenon, and even neon. For the 
diatomic molecules and methane the agreement 
is as expected appreciably less good. The fact 
that oxygen deviates considerably more than 
nitrogen and methane is interesting. This may 
possibly have some connection with the fact that 
the gaseous molecule Oz is in a triplet state. The 
molar entropy consequently contains an elec- 
tronic contribution R log 3, which presumably 
persists in the liquid, but has disappeared‘ at the 
lowest temperatures at which thermal measure- 
ments have been made on crystalline oxygen. 
Whether it disappears by formation of O, mole- 


cules or otherwise is not known, nor in what 


temperature range it disappears, but it may be 
expected to exert a marked effect on the equilib- 
rium properties of the solid phase. 


11. ENTROPY OF FUSION 


The molar heats of fusion L; divided by the gas 
constant R are given in row 20 of Table I and the 
molar entropies of fusion divided by R in row 21. 
The agreement between the values for the 
monatomic molecules is good. 


12. TRIPLE-POINT PRESSURES 


The triple-point pressures P; of the monatomic 
species are given in row 22 of Table I and the 
ratios 100P,/P. in row 23. The values for argon, 
krypton, and xenon, but not neon, agree within 
+4 percent. Since an increase of 4 percent in the 
equilibrium vapor pressure of the crystal, or of 
the liquid, corresponds to an increase in absolute 
temperature of only 0.4 percent, this is a much 
more sensitive test than that of Section 10. 


13. EXPANSION ON FUSION 


In rows 24 and 25 of Table I are given the 
molar volumes v; and », of liquid and solid, 
respectively, at the triple point for the monatomic 
species. In row 26 are given the ratios of v; to 2,. 
As previously pointed out by Clusius and 
Weigand® the uniformity of these ratios is 
remarkable. 


4W. Giauque and H. Johnston, J. Am. Chem. Soc. 51, 
2300 (1929). 

5K. Clusius and K. Weigand, Zeits. f. Physik. Chemie 
B46, 1 (1940). \ 
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14. HEAT CAPACITY OF LIQUID 


According to the principle of corresponding 
states the molar heat capacity at corresponding 
temperatures should be the same for monatomic 
molecules. That this is at least roughly true is 
shown by the values given in row 27 of Table I. 
The application of the principle to diatomic 
molecules is less simple. One of the assumptions 
used to deduce the principle is that the molecules 
should rotate freely in the liquid. If this assump- 
tion is accurately fulfilled, then the rotational 
degrees of freedom should make precisely the 
same contribution to each of the thermodynamic 
functions in the liquid as in the gas. In particular 
the rotational contribution to the molar heat 
capacity of the diatomic species should be just R. 
Hence the total molar heat capacity of the 
diatomic species should exceed that of the 
monatomic species by just R. Actually the excess, 
as shown by the values in row 27 of Table I, is 
greater than R. This is understandable if there is 
a small restriction of the rotation increasing with 
decreasing temperature. 


15. EXTENSION OF PRINCIPLE TO SURFACES 


In order to extend the principle of corre- 
sponding states to surfaces, one needs to revise 
assumption III as follows: ‘‘The intramolecular 
degrees of freedom are assumed to be completely 
independent of the volume per molecule and of 
the distance of the molecule from the surface.” 
From this assumption it follows that the excess 
free energy F, due to an interface of area A con- 
taining NV, molecules is of the form — 


F, kT 


A 


where y, is a complicated but universal function 
of kT/eo and A/N,R,?. It is also evident that at 
equilibrium between liquid and vapor A/N,R¢ 
is like v,/Ro* a universal function of T/T. conse- 
quently F,Ro’/kTA is a universal function of 
T/T.. But for a single component system the 
surface tension y is equal to F,/A. Hence yRo?/kT 
is a universal function of T/T. Finally, using the 
fact that v./Ro® has a universal value, one con- 
cludes that yv.!/kT and so also yv.!/kT, are uni- 
versal functions of 7/7. How well this prediction 
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is borne out by the experimental data will be 
examined in the two succeeding sections. 


16. TEMPERATURE DEPENDENCE OF 
SURFACE TENSION ; 


Since the surface tension of a liquid decreases 
with increasing temperature and vanishes at the 
critical point, the simplest possible form of 
empirical relation between y and T is 


( T. 


where 7¢ is a constant. | shall show that this rela- 
tion with r= 2/9 in fact adequately represents the 
experimental data, but for the moment I leave 
the value of 7 unspecified. Van der Waals sug- 
gested a formula such as (16.1) with r=0.234 and 
Ferguson® has proposed r=0.21+0.015. 

One of the earliest and best known empirical 
relations relating to the temperature dependence 
of y is that of Eétvos, namely 


yui«(1—T/T.). 


(16.1) 


(16.2) 


This relation was however found to be inexact 
and various modifications of it have been pro- 
posed, of which the most satisfactory is that of 


‘Katayama,’ namely 


yy t«(1—T/T.), (16.3) 
where y is defined by 
y =1/v,—1/2., 


We = (pi— Pg) / 


(16.4) 
(16.5) 


Katayama showed that the experimental data 
for various organic substances were in excellent 
accord with his relation (16:3).. How accurately 
this relation fits the data for argon, nitrogen, etc. 
can only be confirmed by trial, but anticipating 


‘such an examination I shall tentatively and 


provisionally assume that the relation (16.3) 
holds for these substances. 

By substitution of (16.5) into formula (6.4), 
one sees that this can be expressed in the form 


y« (1—T/T,)}. (16.6) 


6 A. Ferguson, Trans. Faraday Soc. 19, 407 (1923); Proc. 
‘Phys. Soc. 52, 759 (1940). 

7 Katayama, Science Reports Téhoku Imperial Uni- 
versity 4, 373 (1916). 
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TABLE II. y=yo(1—T/T-) "9. 


Ne A 


N2 Oz 


T.=44.8°K T-=150.7°K 
yo= 15.1 dyne/cm vyo= 36.31 dyne/cm 


¥ 
dyne/cm dyne/cm 
obs. 
13.19 
12.68 
11.91 


calc. 
13.16 
12.67 
11.95 


V.=75.3 cm*/mole 
Yo VATA“ 
=4.3 erg deg.—! mole! 


V.=41.7 cm’/mole 
vo VAT“! 
= 4.05 erg deg.~! mole! 


T.=126.0°K T.=154.3°K 
yo= 28.4 dyne/cm vo= 38.4 dyne/cm 


dyne/cm dyne/cm 
calc. 

18.34 
17.02 
15.72 
14.44 
13.17 


V.=74.5 cm’/mole 
Yo Te! 
=4.4 erg deg.~! mole“! 


obs. 

10.53 
9.39 
8.27 
7.20 
6.16 


calc. 

10.54 

“9.40 
8.29 
7.20 
6.14 


= 90.2 cm?/mole 
Yo VATA 
=4.5 erg deg.—! mole! 


Now eliminating y between (16.3) and (16.6), one 
obtains 


ye (1—7/T), 


(16.7) 


which can also be written 


(16.8) 


in agreement with (16.1) with r=2/9. 

To recapitulate, for argon, nitrogen, and other 
substances under consideration the validity of 
formula (16.6) has already been established in 
Section 6. Taking account of this, formula (16.8) 
and Katayama’s formula (16.3) are mutually 
inter-dependent. Verification of either implies the 
accuracy of the other and conversely. I shall now 
show that formula (16.8) holds satisfactorily for 
the substances under discussion for which experi- 
mental data exist, namely neon, argon, nitrogen, 
and oxygen. ‘ 

In Table II at the top, are given the experi- 
mental value of the critical temperature T, and 
the arbitrarily assumed value of yo for each 
substance. Below this the values of y calculated 
according to (16.8) are compared with the ob- 
served values over a range of temperatures. The 
agreement is good throughout and could hardly 
be improved by amending the value 2/9 assumed 
for r. 

If now instead of eliminating y between (16.3) 
and (16.6) to obtain (16.7), one eliminates T/T. 
one obtains 


s=3(1+7) =33. (16.9) 


This is in definite disagreement with the better 


known relation of McLeod? according to which s 
is 4. Actually McLeod based his relation on half 
a dozen organic compounds, but a mere glance at 
his figures is sufficient to show that y in fact 
varies as some power of y less than 4. 

It is outside the scope of the present paper to 
discuss the best value of s to represent the over-all 
average behavior of organic compounds. I 
merely wish to emphasize that for the substances 
having the simplest and most symmetrical mole- 
cules the value of s is not 4 but 33. Any molecular 
theory or model of the simplest molecules 
leading to s=4 is therefore spurious. 


17. ABSOLUTE VALUE OF SURFACE TENSION 


In Section 15 it was shown that according to 
the principle of corresponding states yv.!7- 
should be a universal function of 7/T7.. Seeing 
that for the substances under consideration the 
dependence of y on the temperature is repre- 
sented by formula (16.8), it follows that yov.t7 
should have a universal value. How nearly this is 
the case is shown at the bottom of Table II. The 
agreement between argon, nitrogen, and oxygen 
is remarkable; that of neon is less good, pre- 
sumably owing to quantal effects. 


18. DEVIATIONS FROM THE PRINCIPLE 


In the column of Table I preceding neon and 
labelled ‘‘Ideal,’”’ are given values of those 
quantities which should be common to all sub- 


< McLeod, Trans. Faraday Soc. 19, 38 (1923). 
® Fowler and Guggenheim, Statistical Thermodynamics 
(Cambridge University Press New York, 1939), p. 451. 
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T T 
°K Ss. 
24.8 5. 70.0 35 
25.7 - 5. 75.0 0 
26.6 80.0 
27.4 475 4.69 85.0 5 
28.3 4.45 4.44 90.0 23 
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stances conforming to the principle of corre- 
sponding states. It has been shown that argon, 
krypton, and xenon conform with a high degree of 
accuracy to the principle in all their properties; 
that nitrogen, oxygen, carbon monoxide, and 
methane also conform with fair accuracy in the 
gaseous and liquid states but not in the solid 
state. The behavior corresponding to argon, 
krypton, and xenon should be considered as an 
ideal one to which other substances may be ex- 
pected to conform more or less according to the 
extent to which the assumptions of Section 2 are 
or are not satisfied. Departures of a particular 
substance from this ideal behavior must not be 
regarded as unsatisfactory flaws in the principle, 
but rather as giving interesting information con- 
cerning fundamental differences between the 
molecules of the substance in question and the 
ideal molecule visualized in Section 2. The varied 
effects of such deviations from the ideal have 
been discussed by others, notably by Hilde- 
brand,'® by Bernal," and by Pitzer.! 


APPENDIX 


The following sources of experimental data 


have been used: 


1. Critical data, densities of coexistent liquid and vapor. 
Neon: Mathias, Crommelin, and Onnes, Comptes rendus 
175, 933 (1922). 
Mathias and Crommelin, Comm. Leiden Supp. 52 
(1924). 
Argon: Onnes and Crommelin, Comm. Leiden No. 131a 
(1912). 
Krypton: Mathias, Crommelin, and Meihuizen, Comptes 
rendus 204, 630 (1937). 
Xenon: Patterson, Cripps, and Gray, Proc. Roy. Soc. 
86, 579 (1912). 
Nitrogen: Mathias, Onnes, and Crommelin, Comm. 
Leiden No. 145c (1914). 
Oxygen: Mathias and Onnes, Comm. Leiden No. 117 
(1911). 
Carbon monoxide: Crommelin, Bijleveld, and Grigg, 
Comm. Leiden 20 No. 221b (1932). 
Methane: Cardoso, J. Chim. Phys. 13, 312 (1915). 


10]. Hildebrand, J. Chem. Phys. 7, 233 (1939). 
1938) Bernal, Address to Chemical Society, February 3, 
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2. Boyle points. 
Neon, argon, nitrogen: Holborn and Otto, Zeits f. Physik 
33, 9 (1925). 
Carbon monoxide: rough interpolation of data of Bartlett, 
J. Am. Chem. Soc. 52, 1374 (1930). 
Methane: Keyes and Burkes, J. Am. Chem. Soc. 49, 1403 
(1924). 


. Virial coefficients. 
Holborn and Otto, Zeits f. Physik 33, 5 (1925). 


. Vapor pressures. 
Landolt-Bérnstein tables. 


. Heats of evaporation and of fusion, triple-points, heat 
capacities of liquids. , 
Neon: Clusius, Zeits. f. Physik. Chemie B4, 1 (1929). 
Clusius, and Riccoboni, Zeits. f. Physik. Chemie 
B38, 81 (1937). 
Argon: Eucken and Hauck, Zeits. f. Physik. Chemie 134, 
161 (1928). 
Clusius, Zeits. f. Physik. Chemie 31, 459 (1936). 
Frank and Clusius, Zeits. f. Physik. Chemie B42, 
395 (1939). 
Krypton: Clusius, Kruis, and Konnertz, Ann. d. Physik 
33, 642 (1938). 
Xenon: Clusius and Riccoboni, Zeits. f. Physik. Chemie 
B38, 81 (1937). 
Nitrogen: Giauque and Clayton, J. Am. Chem. Soc. 55, 
4875 (1933). 
Oxygen: Giauque and Johnston, J. Am. Chem. Soc, 51, 
2300 (1929). 
Carbon monoxide: Clayton and Giauque, J. Am. Chem. 
Soc. 54, 2610 (1932). 
Methane: Frank and Clusius, Zeits. f. Physik. Chemie 
B36, 291 (1937). 


6. Coefficients of thermal expansion. 
Baly and Donnan, J. Chem. Soc. 81, 907 (1902). 


. Densities of liquid and crystal at triple-point. 
Clusius and Weigand, Zeits. f. Physik. Chemie B46 
(1940). 


. Surface tensions. 
Neon: van Urk, Keesom, and Nijhoff, Comm. Leiden 
No. 182b (1926). 
Argon: Baly and Donnan, corrected by Rudorf, Ann. d. 
Physik [4] 29, 764 (1909). 
Nitrogen, oxygen: Baly and Donnan, Trans, Chem. Soc. 
81, 907 (1902). 
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A simple method for the determination of bond energies and internuclear distances is de- 


veloped based on old quantum theory models modified to include the effect of zero-point 
vibrational energy and a type of resonance (orbital pulsations) assumed to be associated with 
nuclear vibrational motion. The virial theorem is employed in evaluating the proposed models. 
The method is applied to the hydrogen molecule and the ‘molecule ion’? (H2*). Excellent 
agreement with experiment is achieved. A possible correlation of the old quantum theory and 


wave mechanics is indicated. 


N the days of the old quantum theory, the 
procedure used in studying atomic and mo- 
lecular systems was to seek out the so-called 
multiply or conditionally periodic motions of the 
particles of the system on the basis of a Fourier 
analysis of some assumed classical model set up 
in conformity with the Rutherford-Bohr concept 

‘of the structure of matter." The models for one- 
electron atoms (H, He’, etc.) and systems which 
could be approximated as such proved eminently 
successful except in some of the finer details, and 
a great deal of the vast accumulation of spectro- 
scopic data was satisfactorily interpreted on this 
basis. On the other hand, the old quantum theory 
failed completely in its treatments of molecules 
and many electron atoms.’ The newer quantum 
mechanics, however, has achieved outstanding 
success in quantitative treatments of the simplest 
systems of this character (He and He-like atoms, 
He, H2*),4-* and in approximate treatments of 
numerous other problems relating to the struc- 
ture and dynamics of matter. As a consequence, 
the old quantum theory is now considered to 
have only historical interest and the structure of 
matter is generally considered as adequately 
described by the quantum mechanical-proba- 
bility distribution functions. 

The theory of wave mechanics, of course, has 

a number of advantages over the older theory, 
one in particular being its straightforward and 


' Rutherford, Phil. Mag. 21, 669 (1911). 

2.N. Bohr, Phil. Mag. 26, 1 (1913). 

3 J.H. VanVleck, Bull. Nat. Research Council 10, (1926). 
(1933) M. James and A. J. Coolidge, J. Chem. Phys. 1, 825 
a eS Hylleraas, Zeits. f. Physik 65, 209 (1930); 71, 

6G, Jaffe, Zeits. f. Physik 87, 535 (1934). 
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unambiguous methods of quantization. It has one 
rather subtle advantage, however, which may 
have contributed more to its success than is 
generally supposed. It is less exacting than the 
older theory in its description of the dynamics of 
matter and the various factors contributing to 
structure. One may, for example, obtain by the 
variation method a very accurate treatment of a 
particular system without the necessity of taking 
into consideration the various contributing phys- 
ical factors, the mathematical quantities em- 
ployed in the wave functions being more or less 
arbitrarily chosen. This is clearly illustrated by 
considering the status of the resonance phe- 
nomenon in the most accurate quantum me- 
chanical treatments, aside from one-electron 
atoms, yet obtained; although most authorities 
agree that the so-called quantum mechanical 
resonance phenomenon contributes greatly to the 
strength of the chemical bond, the treatment by 
James and Coolidge‘ of the Hz molecule makes no 
use of resonance or exchange integrals, or recog- 
nizable terms accounting for this factor. While 
the old quantum theory models could be altered 
through various limits, unless the correct model 
were chosen and all important contributing 
periodic motions taken into account, agreement 
with experimental observations obviously was 
not achieved. 

In this paper, a reconsideration of old quantum 
theory models is undertaken in light of two types 
of motion which were not recognized in the days 
of the old quantum theory, namely (1), resonance 
and (2), ground state nuclear vibrations both 
factors becoming evident from quantum me- 
chanical analyses. In this-enterprise, we recognize 


945 


that present concepts hardly justify a considera- 
tion of classical models in atomic and molecular 
structure. However, if successful, we stand to 
gain materially in simplicity of treatment and in 
clarity of understanding from the structural and 
dynamical viewpoints. Obviously, success in such 
attempts would require some revision in quantum 
mechanics particularly its presently accepted 
philosophical interpretations, although the theory 
as a whole contains far too many points in 
agreement with experiment to be fundamentally 
affected. 


PROPOSED MODEL 


In the old quantum theory method, the working 
model of a particular molecule was set up on the 
assumption that the nuclei would be at rest with 
respect to each other as a result of an equilibrium 
of forces, and that the electron orbitals would 
occupy such stationary positions as to provide the 
necessary stability. We have learned since that 
even in the ground state the nuclei are not at rest 
with respect to each other, but possess a finite 
zero-point energy by no means negligible. Conse- 
quently, if we are to assume that the system in 
its ground state will be determined by the equi- 
librium distance giving a minimum potential or 
total energy (the latter is used in the present 
discussion), the considerations by Slater based on 
the virial theorem’ show as applied to the 
classical models that the nuclear equilibrium dis- 
tance will be greater than that determined by the 
old quantum theory model of rigid nuclei. In 
order that this may be true, however, equilibrium 
must be established not by a balancing of forces 
on the nuclei at every instant of time but by a 
statistical balancing over a period of time. Thus 
at one instant the nuclei will be receding from 
each other as a result of a net repulsion at some 
particular phase, while in the opposite phase the 
reverse will be true. What then will be the effect 
of these nuclear oscillations on the electron 
orbitals? The nuclear vibrations must, of course, 
be brought about by a periodic exchange of mo- 
mentum between the electrons and the nuclei 
occurring in such a manner as to maintain a 
constant total energy and angular momentum. 
Considering the relative masses of the electrons 


7J. C. Slater, J. Chem. Phys. 1, 687 (1933). 
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and the nuclei, we might expect therefore that in 
some cases relatively large electron orbital pulsa- 
tions or oscillations will take place during a cycle 
of nuclear vibration. Thus a type of resonance is 
indicated which is intimately associated with the 
vibrational energy. While this type of resonance 
differs somewhat from the quantum mechanical 
resonance phenomenon, we may discover that the 
two types have much in common. 

As suggested by London,*® many chemical re- 
actions are “‘adiabatic’’ in the sense that they 
involve no electronic transitions. In such cases, 
the character of the orbit, i.e., whether circular or 
elliptical in the view of the old quantum theory, 
will be relatively unchanged in the passage from 
the atom to the molecule. This probably will be 
nearer to the true state of affairs in simple-mole- 
cules than in complex ones, for in the latter type 
the complex spacial distribution of charge may 
cause considerable orbital change or distortion 
despite the absence of monochromatic radiation 
in the formation of the chemical bond. Such 
factors, however, need not be considered in the 
present study since we here deal only with the 
simplest molecules (Hz and H2*). Thus we shall 
assume in discussing the ground state of these 
simple molecules that the electron orbitals are 
essentially though not exactly circular. Since, 
however, the frequency of nuclear vibration will 
be small relative to the electron orbital fre- 
quencies, we shall assume circular orbits for 
methods of calculation, and neglect the small 
energies associated with the frequency of pulsa- 
tion as well as electron spin and relativity effects. 
Indirectly, of course, we shall expect that the 
assumed orbital pulsations caused by nuclear 
vibration will contribute largely to the binding 
energy not because the pulsation frequencies 
themselves are high but simply because the 
pulsating orbits will allow stabilization of the 
nuclei at such separations as to lead to a greater 
binding energy, for example by an increase in the 
ratio | V| to T (V-potential energy, and 7-kinetic 
energy).” 


METHOD OF ANALYSIS 


The virial theorem applied as mentioned above 
with a minimum of refinements, i.e., by intro- 


8 F. London, Zeits. f. Elecktrochemie 35, 552 (1929). 
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ducing only the electronic principle quantum 
number mz, provides a surprisingly accurate 
method for investigating the proposed model, 
although quite naturally this simple treatment 
does not claim the accuracy inherent in the most 
refined quantum mechanical treatments of the 
systems considered here.*~* As shown by Slater,’ 
the total energy E for a given interatomic dis- 
tance R will be equal tothe Energy E’ determined 
‘on the assumption that the system is under no 
constraint (i.e., by the virial theorem relation 
between T and V), plus the virial of external 
constraint. Hence by Slater’s method for a 
diatomic molecule 


E=T+V=E'—(RF), (1) 


F, the external force of constraint required for 
equilibrium, being positive for attraction and 
negative for repulsion. While we shall adopt 
Slater’s procedure of regarding the nuclei as fixed 
at a given value of K “y -> external constraint 
sufficient for equilibri will be necessary in 
the model under consic on in order to obtain 
the true virial of constraint to call to mind at 
various times our premise that the nuclei actually 
are equilibrated by a pulsation or oscillation of 
orbitals. This consideration would be unnecessary 
except for the fact that the orbital dimensions 
and positions will depend to some extent on the 
particular phase of motion involved. On the 
average, we have for a diatomic molecule under 
inverse square forces 


[= —E'—}(RF)=—E-(RP), 
=2E’ =2E+(RP).- (2) 


In obtaining (RF) and E’, however, the average 

orbital radii must be found by averaging over a 

complete cycle of motion, or more simply over 

the two opposite phases at precisely the equi- 
librium position. Hence, in the phase of inward 
motion at the equilibrium position 


Ti = 4RF,, 
Y= 2Ey’, 
while in the opposite phase at this position 


T2.= — —E;'+3RFi, 


since at this position equal but oppositely directed 
forces of constraint would have to be applied to 


(3a) 
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equilibrate the nuclei. The conditions of (2) and 
(3) will be fulfilled only for particular average 
values of the orbital radii. In other words, 
different energies actually are obtained by (3a) 
and (3b), but this is due to the neglect, in each 
separate calculation of the virial of constraint, of 
the actual motions of the nuclei which are moving 
at maximum velocity at the equilibrium position 
but in opposite directions in opposite phases of 
motion. Since the correction factors in the two 
phases at the equilibrium distance are equal in 
magnitude and opposite in sign, clearly a simple 
arithmetical mean average of the values E; and 
E; calculated from Eqs. (3a) and (3b), respec- 
tively, cancels the correction factors and gives the 
correct energy and the virial of constraint (a 
time-average value). 


THE HYDROGEN MOLECULE ION (H:2*) 


First we recognize according to the assumed 
model that the electron orbital at any instant 
should occupy a position such as to apply an 
equal force on each proton. Further, if the orbit 
is to be circular, it must be so oriented with 
respect to the protons that the resultant force on 
the electron is always directed along the radius 
of the orbit. These conditions are fulfilled by the 


_ model in which the orbital lies in a plane perpen- 


dicular to a line joining the protons, and either 
midway between them or at a position on either 
side of the midpoint such that 


(e?/R1?) sin (e?/R2?) sin 


where e is the charge on the electron, R; is the 
distance from the electron to proton 1, Rez is the 
distance from the electron to proton 2, and 6 and 
¢ are the corresponding angles between the orbital 
plane and the lines R; and Ro, respectively. The 
solution as expected from symmetry considera- 
tions corresponds to 6=¢. This is similar to the 
so-called circular model studied in the days of the 
old quantum theory.* ® 

For the symmetrical model, 0=¢, Eq. (1) gives 


E=(-—e/2r)(f+g), (4) 
where 
f=2 cos @—1/2 tan 6, 
g=(R/r)(sin cos? 6—1/4 tan? @), 


9N. Bohr, Phil Mag. 26, 857 (1913). 
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TABLE I, Values of E for various assumed values of R/r for 
(e?/2a9= 13.538 ev). 


—E (ev) 


15.81 
15.94 
16.03 
16.10 
16.13 
16.149 
16.154 
16.159 
16.161 
16.148 
16.11 
16.03 
15.86 


g 


47° 04’ 1.1625 
47° 43’ 1.1731 
48° 22’ 1.1820 
48° 59’ 1.1902 
49° 36’ 1.1968 
50° 12’ 
50° 39 6. 

50° 46’ 1.2078 
50° 55’ 1.2089 
51° 00’ 1.2091 
51° 20’ 1.2111 
52° 26’ 1.2160 
53° 28’ 1.2180 


ann wm 


~ 
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and r is the radius of the orbit. Employing the 
quantum condition v=nh/2xmr (v=electron or- 
bital velocity, m= mass of the electron), Eq. (3a) 
gives 71 = d0/f(a9=.5285A — the Bohr radius), and 
(3b) gives r2=ao/(f+g). In calculating the 
energy from (4) by use of a single value of 7, a 
root mean square average of 7; and 72 should be 
taken since the energy turns out to be propor- 
tional to 1/r?, f and g being inversely proportional 
to r. Taking therefore 


(5) 


r= 


Eq. (4) becomes 


— 2 
2a0 [(f+g)?+f?]}! 


Table I lists values of E for various values of the 
ratio R/r;a minimum value of — 16.16 ev occurs 
at R/r=2.46, r=.535A and R=1.32A. This cor- 
responds to a bond energy Dp» of 2.62. ev; the 
experimental value is 2.639 +.002 ev (our method 
of calculation, of course, includes the zero-point 
vibrational energy). 

The discrepancy in R is only apparent; the 
value 1.32A actually corresponds to one ex- 
tremity of nuclear motion. The reason for this is 
that we used the same value of R/r in obtaining 
both 7; and 72 which condition applies only at one 
extremity of motion where the nuclei are mo- 
mentarily at rest. We may with equal simplicity 
calculate both the true equilibrium position and 
the opposite extremity of nuclear vibration by 
careful reference to the assumed model. In these 
determinations, however, it is necessary to em- 
ploy the value of E determined by (6) because 
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apparently no minimizing conditions appear in 
the method employed for these other positions of 
nuclear vibration. The reason for this will become 
evident as we explore the various possible separa- 
tions of the nuclei by reference to the model in 
use. 

If we simply replace the term f+g by f—g in 
(6), a value of R=.73A is obtained at R/r=1.36 
for constant E (— 16.16 ev). This conforms to our 
model for the opposite extremity of nuclear 
vibration remembering that in this position the 
nuclei are in a state of repulsion corresponding to 
a negative virial of constraint. 

Now at the equilibrium position we are dealing 
with two opposite phases of motion so that the 
value R/r will be different for each phase. By 


' assuming a value of R instead of R/r, we may 


readily evaluate r; and r2 by successive approxi- 
mations each value of R/r corresponding to a 
given value of f and g, since R/r=2 tan 6. Re- 
quiring again that E=constant (— 16.16 ev), we 
obtain R=1.07A at r1;=.5708A and r2=.4380A. 
The experimental equilibrium distance is 
R=1.06A. If an attempt is made to minimize E 
by reference only to the equilibrium position, the 
value ~ —17.0 ev is obtained at R#1.3A. How- 
ever, a displacement of the nuclei by a small 
amount on either side of the position R#1.3A 
would then require an increase in E. Since E isa 
constant of motion, this would mean that the 
amplitude of nuclear vibration would be zero, a 
condition not in conformity with the assumed 
model. This explains why it is necessary to 
minimize E in the present method for one ex- 
tremity of nuclear vibration. 


THE HYDROGEN MOLECULE 


The circular model proposed by Bohr was also 
applied to Hz in old quantum theory studies. In 
this case, the electrons were considered to occupy 
a common circular orbit at diametrically opposite 
positions, the orbit lying in a plane perpendicular 
to the line joining the nuclei and midway between 
them. Applying the methods of this study to a 
model of this character modified to include 
nuclear vibrations and orbital pulsations, we 
obtain an energy that is too low (too large nega- 
tive energy). One might be led to expect, how- 
ever, that such a system in which the electrons 
are required to move in a common orbit would be 
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Fic. 1. Model of He for a particular phase of motion. 
Orbitals are considered to pulsate (r to oscillate), and to 
oscillate between the nuclei A and B. Thus at an opposite 
phase of nuclear motion the positions of electrons 1 and 2 
may be reversed and at the nuclear equilibrium distance the 
orbitals may coincide. 


dynamically unstable; this was one of the argu- 
ments against the circular model in the days of 
the old quantum theory.® 

Let us investigate one possible way by which 
the Hz molecule might be formed from two H 
atoms. We suppose that at some relatively large 
internuclear distance R the two then independent 
systems will exert a force of polarization on each 
other such as to cause the orbitals to line up 
parallel to each other in planes perpendicular to 
the line AB (Fig. 1) joining the nuclei. The 
energy for a given value of R then would be a 
minimum when the electrons both rotating in the 
same sense were always in the same plane through 
the line AB (the plane rotating with the elec- 
trons) but at opposite sides of this line. Due to 
the identity of the protons, and of the electrons, 
we assume that the system will be completely 
symmetrical as long as the two atoms remain 
under the influence of each other, i.e., each 
orbital will occupy the same position relative to 
its own proton and that of the other orbital, and 
the orbital radii will be the same at any instant. 
Under this configuration there will exist a net 
attraction so that the atoms, in which the orbitals 


will be drawn away from their respective nuclei 
toward the center of the system, will approach 
along the line AB until the limits of the molecule 
are reached. At this stage the nuclear oscillations 
and electron orbital pulsations presumably will 
set in if not earlier. The electron orbitals may 
possibly oscillate back and forth between the 
nuclei, and in case they approach close enough to 
the center position of AB, the orbitals may cross 
over at the center point (@=@), and exchange 
positions periodically during each cycle or half- 
cycle of nuclear vibration. 

For the above model we may again apply an 
equation for E similar to (4) neglecting as before 
spin and relativity effects, and energies associated 
with frequencies of orbital pulsations and oscilla- 
tions. The definitions of f and g in (4) now become 


COs @ 1 
f =2 cos 6+2 cos 
2 tan 0+tan 


and 


R 
g -—(sin 6 cos? 6+sin ¢ cos? 
r 


1 
@ 
(tan 6+tan ¢)? 


The definitions of @ and ¢ for each electron are the 
same as in H»* except that they are not necessarily 
equal in Ho, a is the angle between a given orbital 
and a line passing through the two electrons, and 
r is the radius of either orbit (Fig. 1). From the 
quantum conditions, and (3a) and (3b), we obtain 
71 =2a0/f and r2=2ao/(f+g). Using again a root 
mean square average for 7, (1) becomes 


(2)'(f+g)?f 
4ay 


differing in form from (4) only by a factor 2 in 
the denominator. 

Table II gives values of the total energy and 
other factors employed in the calculations for 
various assumed values of R/r. From these data a 
minimum total energy of —31.55 ev is seen to 
occur at R/r=1.24, r=.495A and R=.61A. The 
value of Dy obtained by this method is thus 
4.47 ev. The experimental value is given as 
4.454+.013 ev. 


(8) 
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As in H,* the value of R given above corre- 
sponds to an extremity of nuclear vibration which 
surprisingly in this case appears to be the inward 
rather than the outward extremity. That this is 
the case is readily verified when we find that even 
in the special case 6=¢ and a=0, the energy 
increases (becomes less negative) for all smaller 
values of R. However, at R=.61A the model 
6=¢ gives almost exactly the same energy value. 
Hence we see immediately that the possibility 
exists of orbital interchange mentioned above. 
While the minimum energy for the case 6¥#@ 
comes at R=0.61A, we find on investigating the 
model @=¢ (by taking simply r= 2ao/f) that the 
minimum value (£:2—35.3 ev) occurs at 
R/r=1.6, or R&.85A. It seems logical to assume 
therefore that upon reaching the position R= .61A 
in the formation of the molecule, orbital pulsa- 
tions set in by the electrons first moving toward 
a common orbital, the orbitals suddenly ex- 
panding as R tends to decrease below .61A thus 
allowing a nuclear repulsion to set in. The 
tendency then is for the nuclear distance to 
increase to R~.85A. The tendency for a further 
decrease in energy below — 31.55 ev is, however, 
upset by the passage of the orbits over each other 
to assume again the condition 6#¢. The repulsion 
toward the outward position probably corre- 
sponds, however, to the model @=¢ where 
re=2a0/(f—g) and r1=2a/f in place of the 
values used in deriving (8). At least we might 
expect that at the exact equilibrium position the 
model 6=¢ will hold. Let us, therefore, seek a 
solution of the equilibrium position in a manner 
similar to that employed in He", i.e., by taking E 


TABLE IT. Variation of E with R/r for He for the condition 
that Using r1;=2a0/f and r2=2a0/(f+g). 


R/r 0 g —E (ev) 


1.154 30° 30° — — 29.79 
1.20 30°45’ 31°13’ 21") 2.0959 .0748 31.33 
1.22 26°30’ 35°48’ 6°22’ 2.0954 .0856 31.54 
1.23. 25°15’ 37°10’ = =8° 15’ 2.0945 .0870 31.55 
1.24 24°03’ 38°26’ 9°51’ 2.0937 .0881 31.55* 
1.25 23°10’ 39°25’ 11° 9’ 2.0931 .0892 31.55 
1.28 20°52’ 41°57’ 14°30’ 2.0910 .0899 31.52 
1.30 19°40’ 43°18’ 16°18’ 2.0897 .0910 31.50 
1.34 17°45’ 45°34’ 19°17’ 2.0868 .0934 31.46 
140 15°27’ 48°20’ 22°57’ 2.0837 .0945 31.40 
145 13°55’ 50°15’ 25°30’ 2.0793 .0948 31.27 
1.50 12°39’ 51°54’ 27°44’ 2.0762 .0956 31.20 
1.60 10°35’ 54°42’ 31°31’ 2.0700 .0946 30.98 


* At R/r =1.24, r =.4948A and R=.61A, 
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constant (—31.55 ev) and average the energies 
E,; and £, obtained from (3) by first taking 
r=2ao/f and then r=2a9/(f—g) for a particular 
value of R. Table III gives the results of this 
calculation from which we see that — E= 31.55 ev 
at approximately .72A and .88A for R. 


TABLE III. Values of E vs. R for He for the condition 6=4¢, 
a=0. Using r2=2a0/(f—g) and r1;=2a0/f. 


—E (ev) 
(ev) -—E:(ev) =—E:—E:/2 


R (A) v1 (A) ro (A) 
.70 .5528 5111 29.33 33.46 31.40 
72 .5610 29.19 33.93 31.56 
74 .5693 .5150 29.02 34.16 31.59 
.80 .5934 5213 28.54 34.94 31.69 
84 .6095 .5262 28.21 35.35 31.78 
85 .6127 .5285 28.12 35.30 31.71 
.87 .6192 .5313 27.95 35.25 31.60 
.90 .6340 .5365 27.65 35.15 31.40 


—E =31.55 ev at R&.72 and .88A. 


We thus conclude from this method of treat- 
ment that the nuclei in He oscillate between 
R=.61A and R=.88A with an equilibrium dis- 
tance of about .72A. The experimental equi- 
librium distance is given as .74A. 

The determination particularly of the outward 
extremity of nuclear vibration is somewhat arbi- 
trary. Possibly the outward extremity should be 
calculated by averaging between the conditions 
corresponding to 6=¢ and @#¢@. However the 
result cannot be far different since the energy 
thus calculated will be found to increase rapidly 
above R=.88A. A more satisfactory and certainly 
more elaborate treatment would be to solve the 
equations of motion by evaluating gpdg. This 
procedure, however, may prove quite difficult. 
In such a treatment, the results of the virial 
treatment might serve as a useful guide in 
establishing the types of motion involved. 


DISCUSSION 


The results of the present treatments of Hz and 
H,* appear to provide further substantial evi- 
dence in favor of the Rutherford-Bohr models as 
applied to atomic and molecular structure. It is 
not surprising, however, that a classical counter- 
part of wave mechanical treatments should be 
found for molecules as well as for one-electron 
atoms. Also quantum mechanics has drawn 
heavily from classical concepts in almost all 
phases of application, a formal similarity of 
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course being recognized in the formulation of 
wave mechanical operators. By way of providing 
some philosophical justification for our return to 
old quantum theory concepts, the following very 
_ brief discussion is presented suggesting a possible 
basis of correlation between wave mechanics and 
old quantum theory. 

The original proposal by de Broglie,!° which led 
to the formulation by Schrédinger" of the theory 
of wave mechanics, described matter as dual in 
nature exhibiting properties of both waves and 
discréte particles, the particle and wave charac- 
teristics being related through the equation 


(9) 


dX being the wave-length and mv the classical 
momentum. It included also the suggestion that 
stationary states of matter may correspond to 
standing ‘‘matter waves.’ Thus in closed orbits 


the total perimeter of the classical orbit corre-— 


sponds to an integral number of wave-lengths of 
the ‘‘matter wave’ so that we have 


(10) 


=n(n=1, 2, -), 


where ds is an element of length of the classical 
orbit. Substitution of (9) in (10) gives directly the 
Wilson-Sommerfeld quantization rules,” 


where » and g are the generalized momenta and 
cordinates. 

Is it not possible that the observed dual nature 
of matter may be due to a direct association of 


(11) 


10L. de Broglie, Thesis, Paris 1924, Ann. de physique 
[10] 3, 22 (1925). 

1 E. Schrédinger, Ann de physique 79, 316, 489, 734; 80, 
437; 81, 169 (1936). 

22 Wilson, Phil. Mag. 29, 795 (1916). 

18 Sommerfeld, Ann. de Physique 51, 1 (1916). 
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electromagnetic waves and discrete particles both 
being understood in the classical sense? In other 
words in the absorption of an electromagnetic 
wave by matter, the wave may not loose its 
identity but merely exist as an ‘‘atmosphere”’ 
surrounding the particle the latter with its 
electrical charge being forced to “‘ride the wave.’’4 
(We do not doubt that electromagnetic waves are 
substance or at least are associated with sub- 
stance other than gross matter.) Om this basis, 
wave mechanics ‘might be considered to be a 
treatment of the electromagnetic ‘‘atmosphere”’ 
the (particle) distribution of which would then be 
described by the probability distribution func- 
tions of that theory. The classical orbit of the par- 
ticle would be determined by Eqs. (9), (10), and 
(11), i.e., by the forces of association of the parti- 
cle and its electromagnetic wave, the latter exist- 
ing as a standing wave as determined by the par- 
ticular boundary conditions of the system. The 
inner quantum numbers of the old quantum 
theory and wave mechanics as is well known do 
not agree. Possibly the two sets of quantum 
numbers are simply orthogonal in the same sense 
that the electromagnetic field exists in relation to 
an electrical current. (We make no attempt here 
to explain the diamagnetic properties of the 
hydrogen molecule.) On a basis such as this, the 
Wilson-Sommerfeld quantization rules may again 
provide a sound method for quantization at least 
for particles moving in closed orbits; for open 
orbits and closed ones for which the various 
periodic motions may not be readily recognized 
the more elaborate and fundamental methods of 
wave mechanics would be required to provide the 
necessary quantization. 

In conclusion, the author wishes to express his 
thanks to Professor Henry Eyring, Princeton 
University, for encouragement and valuable 
criticism and advice. 


4B. Leibowitz, Phys. Rev. 66, 343 (1944). 
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INTRODUCTION 


HE development of transmutation tech- 
niques has facilitated the study of elements 
hitherto undiscovered, or available only in limited 
quantity in their stable form, or even non- 
existent in a stable form in nature. Such has been 
the case of Element 43, discovery of which was 
claimed! in 1925, but concerning whose properties 
very little was known until the investigations of 
Perrier and Segré®* established conclusively that 
the element did exist in radioactive forms and 
demonstrated, by means of tracer reactions, 
many of its chemical properties. 

Some oxidation-reduction reactions of the 
element were reported by these workers,’ but 
quantitative data regarding its electrochemistry 
have been lacking. In this paper are reported the 
results of experiments on the electrodeposition 
of Element 43 from aqueous solutions, and a de- 
termination of the standard potential of the ele- 
ment between its metallic and highest valence 
state has been made. Some other miscellaneous 
observations relating to the electrochemical prop- 
erties of the element are also recorded. 


PRINCIPLE OF THE METHOD 


The method used to determine the standard 
electrode potential of Element 43 was one de- 
veloped by Hevesy and Paneth,‘ and employed 
by them and others® * for determining deposition 
potentials of the natural radio-elements, espe- 
cially those obtainable in limited quantity. Based 
on the assumption that the Nernst law is valid* 
at extremely low activities, the method enables 


* Present address: Department of Chemistry, Uni- 
versity of Minnesota, Minneapolis, Minnesota. 

1W. Noddack, I. Tacke, and O. Berg, Sitz. Preuss. 
Akad. Wiss. 19, 400 (1925). 

2 C. Perrier and E. Segré, J. Chem. Phys. 5, 712 (1937). 

3C. Perrier and E. Segré, J. Chem. Phys. 7, 155 (1939). 

4G. Hevesy and F. Paneth, Wien. Ber. 122, 1049 (1913); 
123, 1619 (1914). See also Physik. Zeits. 15, 801 (1914). 

5 F, Joliot, J. chim. phys. 27, 119 (1930). 

5 M. Haissinsky, J. chim. phys. 32, 116 (1935). 

* Haissinsky (reference 6) showed that for bismuth the 
law was valid down to 10-” molar Bit***. 


the deposition potential of the element in ques- 
tion to be determined with fair accuracy. This is 
done by measuring the quantity of radio-element 
deposited on an electrode at various electrode 
potentials, using the radioactivity of the element 
as a means of estimating quantities either on the 
electrodes or in solution. A plot, then, of the 
amount deposited in unit time vs. the electrode 
potential, will indicate the potential at which 
deposition of the element begins. Assuming that 
this potential is the true reversible potential for 
the electrode reaction in question, then if ac- 
tivities are known, the Nernst equation may be 
used to calculate E°. 

In the case of Element 43, it is reported* that 
from an acidic solution, deposition occurs on the 
cathode, presumably with the liberation of the 
free metal. The fact of the cathodic deposition 
has been verified, and we have also shown con- 
clusively* that no anodic deposition occurs from 
acidic solutions. 

The method used, then, has been to measure 
rates of cathodic deposition of Element 43 from 
weakly acid solutions, using radioactive isotopes 
of the element to follow its course, and from this 
information to calculate the value of the standard 
electrode potential corresponding to the probable 
electro-reduction process. 


MATERIALS AND APPARATUS 
Element 43 


This was produced by bombardment of molyb- 
denum plates, assay 99.6+ percent Mo, with 
14-Mev deuterons in the Massachusetts Institute 
of Technology cyclotron for 30-50 vamp. hours. 

As large amounts of radiomolybdenum were 
also produced in the bombardment, precaution 
was taken in the separation procedure to avoid 


* Conservation experiments showed that, after elec- 
trolysis, the amount of Element 43 left in solution, plus 
what could be dissolved from the cathode, equalled the 
total amount originally placed in the electrolysis cell. 
Tests of the anode showed that no activity could be dis- 
solved from it, and that it was not radioactive. 
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contamination of the active 43 with active 
molybdenum. The possibility that other radio- 
active elements, such as Zr or Cb, were produced 
in the molybdenum plate and might be mistaken 
for Element 43 was considered carefully. From 
the work of Perrier and Segré? as well as from our 
own experiments, it seemed very unlikely that the 
activity, apart from that due to molybdenum, 
could be traced to any element except 43. Ex- 
periments made with the 43 activity on chroma- 
tographic adsorption columns seemed to offer 
additional evidence that this was so. 

The procedure for separating Element 43 
and molybdenum was as follows: The plates, 
weighing from 0.5 to 2 grams, were dissolved in 
approximately 30 ml of aqua regia, excess HNO; 
was added, and the solution evaporated to dry- 
ness under an infra-red lamp. The evaporation 
with HNO; was repeated until the molybdic 
anhydride was perfectly white. It was then taken 
up with a minimum amount of 15 NH,OH, and 
the solution diluted to approximately 150 ml. 
The solution was then acidified with HNO; to 
methyl orange, 1 gram ammonium acetate added, 
and the solution heated to boiling. Next 1N 
AgNO; was added slowly, with stirring, until 
about a 10-percent excess was present above the 
amount required to precipitate the AgeMoQ,. 
The precipitated Ag2MoQ, was filtered and 
washed with 0.5-percent AgNO;. Meantime, 
some inactive (NH4)2MoO, was prepared from a 
portion of inactive molybdenum plate by the 
process just outlined. This was added to the com- 
bined filtrate and washings from the first MoO, 
precipitate until the molybdate ion was in excess. 
Additional AgNO; was then added, and the 
precipitate washed and filtered as before. Excess 
silver ion was removed from the filtrate by 


adding a slight excess of NH,I and filtering off - 


the Agl. (The precipitate may be washed with 
0.2-percent NH,lI solution.) To the filtrate and 
washings was added aqua regia, and the solution 
evaporated. As the volume decreased, HNO; was 
added from time to time. Finally, about 1 ml of 
conc. H2SO, was added and the solution evapo- 
rated to fumes of SO;. Following careful dilution 
with water, the solution was filtered and was then 
ready for use in the electrolysis experiments. The 
reasonable assumption that the element will be in 
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its highest valence state at this point is supported 
by later evidence. 

This procedure, while time-consuming, has 
advantages in that (1) complete removal of 
active molybdenum is assured* (2) no carrier is 
necessary (although no claim is made that 100 
percent of the 43 is recovered) and (3) the ac- 
tivity is obtained in clear solution free of any 
substances which will interfere in future ex- 
periments. 

All samples of 43 used in subsequent experi- 
ments were prepared in the manner described, 
using reagents of the highest purity obtainable, 
and making a double precipitation of silver 
molybdate. 

The electrolysis circuit used for making the 
depositions is shown in Fig. 1. It consisted of a 
means of applying a given potential to a cathode 
upon which the activity was to be deposited. The 
e.m.f. was supplied by a'5.5-volt storage battery 
and was varied by means of a 300-ohm rheostat 
potential divider. The electrolysis cell consisted 
of a test tube 1.5 cm by 10 cm. The electrodes 
were smooth platinum of the Pregl micro type. 
The gauze cathode had roughly the dimensions 
1 cm by 3 cm. The reference electrode was a 
saturated calomel half-cell and was connected to 
the electrolysis vessel by a salt bridge of saturated 
KCI. The geometry of the electrodes and of the 
calomel reference was rigidly maintained by 
wooden blocks. The potentials were determined 
by means of a Leeds and Northrup student-type 
potentiometer. No attempt was made to control 
the temperature at which the runs were made; it 
was generally in the region of 23°C. 


EXPERIMENTAL PROCEDURE 
Measurement of the Activity 


The activity was measured with a Geiger- 
Miiller counter of the liquid type described by 


* Since the separation of radiomolybdenum from 43 is 
very critical, the procedure was checked most carefully. A 
sample of bombarded molybdenum was carried through the 
complete separation without performing any re-precipi- 
tation. The active solution was then divided into two equal 


portions. One portion was ke 7 as a control. Inactive 


ammonium molybdate was added to the second portion, 
and this portion was carried through the entire procedure 
once again. The total activity ms the second portion, after 
making the separation, was only 5 percent less than the 
total activity of the control portion. This could be attri- 
buted to normal sources of error, such as adsorption on the 
various precipitates, in the eee separation. 
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Bale.? The counter was calibrated using a solu- 
tion of uranium nitrate, and was found to have 
an efficiency of about 1 percent. Knowing this, 
and the apparent half-life of 15-16 days of the 
active 43, an approximate value of the 43 concen- 
tration in the runs could be calculated. 

In every case two ml of the solution under 
count was placed in the counter cup using the 
same pipette. All measurements were made on 
solutions of very nearly the same strength, and a 
sufficient number of total counts was taken in 
each case to keep the probable error to about 1 
percent. After each run the counter tube and cup 
were cleaned with dichromate cleaning solution 
and thoroughly rinsed. Background counts were 
taken frequently. 


Preparation of the Electrodes 


The electrodes were first freed of activity by 
contact with warm 6N nitric acid for 15 minutes. 
They were next treated with dichromate cleaning 
solution at 150°C, followed by thorough washing. 
The electrodes were finally dried under an infra- 
red lamp (200°C) and handled only with nickel 
crucible tongs. 


Preparation of the Cell 


The cell was cleaned with hot dichromate 
cleaning solution after each electrolysis. After 
being rinsed it was flamed dry. 


Making the Run 


Five ml of active 43 in sulfuric acid was 
pipetted into the electrolysis cell. The electrodes 
and the calomel reference were then inserted and 
the whole assembly seated into wooden blocks to 
establish a fixed geometry. The potentiometer 
was set at a given fixed value and the cell allowed 
to stand for 15 minutes to make certain that the 
active solution had attained equilibrium with the 
salt bridge. The external e.m.f. was then applied 
and controlled manually by means of the po- 
tential divider so as to maintain the cathode 
potential constant. It was not difficult to control 
this potential manually after the first 15 seconds 
of the electrolysis in the absence of external 
stirring. After the electrolysis had proceeded for 
a definite length of time, a small siphon was 


7W. F. Bale, Radiology 35, 184 (1940). 
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. Fic. 1, Apparatus. 


inserted and the entire contents of the cell drawn 
into a.clean, dry, test tube and thoroughly mixed. 

A count was then taken on the original solution 
and on the solution after electrolysis. This made 
an accurate knowledge of the half life unnecessary. 
The quantity of material deposited on the cathode 
was determined by difference, on the assumption 
that no activity is lost to the anode or to the 
walls of the electrolysis vessel. This assumption 
was verified on five different occasions by a com- 
plete recovery of activity from the cathode. 

The cell and electrodes were then cleaned as 
above, and the same procedure repeated for each 
of the points to be found in the data presented in 
the next section. 


EXPERIMENTAL DATA 


Rate studies were made at four constant 
cathode potentials using portions of a single 
original solution of active 43 in sulfuric acid. The 
solution had a pH of 2.36; this value was sub- 
stantially unchanged in every case after the 
electrolysis. 

The experimental data is summarized in 
Table I. The total time of electrolysis is given in 
column (a). Column (d) lists the net count per 
minute of the original solution as determined by 
the Geiger-Miiller counter, and column (f) lists 
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TABLE I, 


(b) (c) (d) 


Net 


Initial Background (orig.) 


(e) (h) (i) 


Prob- 
able 
error 


% de- 


Final posited 


0.0 volt vs. saturated calomel 


7776 


32 


243 228.1 


214.1 


202.7 


12.4 214.1 


230.5 


194.6 


0.4 volt vs. saturated calomel 


12.8 131.2 


3520 


31 


113.6 


net count per minute of the electrolyzed solution. 
* The difference between the net original and the 
net final activities is the quantity of activity 
deposited on the cathode and is listed in column 
(g). By expressing the cathode deposit as a 
percentage of the original, as in column (h), the 
activities are converted to a unit basis. The 
numerical ratios listed in columns (b), (c), and (e) 
are the ratios of total count to total time that 
were used in determining the activity in counts 
per minute. The probable error in the measure- 


ments is listed in column (i) for the least accurate 
determination in each series. 


RESULTS 


In Fig. 2 are plotted the data of Table I 
showing the change of cathode activity with time 
for each of the four constant potentials at which 
the solutions were electrolyzed. The curves are of 
similar shape to those obtained by both Joliot® 
and Haissinsky® in their studies on the electrode- 
position of polonium. «Mathematical analysis 
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TABLE I.—Continued. 


(c) (d) 


Net 


Background (orig.) 


(e) 


Final (final) ‘Diff. posited -—error 


1705 
130 


13.8 198.0 


13.1 180.9 


216 


90 216 13.3 


0.6 volts vs. saturated calomel 


0.8 volt vs. saturated calomel 


7176 
> 121.8 108.7 72.2 39.9 


62.7 


93 


7628 


* Quantitative recovery of activity from cathode. 


shows that our curves do not follow a definite 
order, and no attempt will be made to interpret 
their general form in terms of the processes 
occurring in solution. However, that diffusion of 
the active ion is the rate determining step in the 
electrolysis process seems probable from similar 
work by Joliot,> and from consideration of the 
somewhat analogous case of electrolysis with a 
dropping mercury electrode. 

The strength of the active cathode deposit is 
directly proportional to the current density* and, 
as a result, Fig. 2 actually represents the change 
in current density with time at various constant 
voltages. If time is now held constant, it is 
possible to establish the relationship between the 
current density or the percentage deposit and the 
voltage of the cathode. This is essentially the 
treatment used in determining decomposition 
voltages, and such a plot would be expected to 


os The current carried by the ions of Element 43 is, of 
course, immeasurably small at the extreme dilution used. 


yield the same type of curve as results from the 
increase in current at the onset of decomposition 
in experiments of that type. 

In Fig. 3 the strength of the cathode deposit is 
plotted against potential for the three different 
time intervals designated by the broken lines in 
Fig. 2.* 

The result, as we would have predicted, is a 
series of straight lines differing in slope but 
arising from a common origin. The common 
origin is the value of the cathode potential at 
which the activity just starts to deposit and 
gives us directly the electrode potential of the 
reaction which the active 43 undergoes in 
solution. 

The best value for the electrode potential is 
taken as —0.1 volt vs. the saturated calomel 
electrode. The negative value refers to the 
polarity of the cathode with respect to the 


* The selection of the time interval has relatively small 
effect on the extrapolated value of E. 
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Fic. 2. Rate curves at constant potential. 


saturated calomel reference. To check this value, 
a 45 minute electrolysis was conducted at —0.2 
volt. The activity was recovered completely 
' from the solution within the experimental accu- 
racy of the measurement. 

The potential of the calomel reference was 
checked against a quinhydrone electrode and the 
value obtained was in good agreement with the 
accepted value of —0.2464 volts. Furthermore, 
this experiment confirmed the polarities assigned 
to each electrode. 

Assuming that the Element 43 is present as an 
anion, the over-all process occurring in solution 
might be written 


It has been shown that the best value for the 
electrode potential of this reaction in a solution 
with a hydrogen ion activity of 10-?:* is —0.10 
volt vs. the saturated calomel ‘electrode. This 
corresponds to a potential of —0.146 volt on the 
normal hydrogen electrode scale. 

Consider the equation further. The activity 
(effective concentration, not radioactivity) of the 
metallic masurium is unity by definition, and the 


8F. MacDougall, Physical agers (The Macmillan 
Company, New York, 1943), p. 


0.146 =E°— 
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activity of the water is essentially unity in these 
dilute solutions. If a reasonably accurate value 
can now be assigned to the activity of the 
permasurate ion, it is possible to apply the 
Nernst equation and calculate the normal elec- 
trode potential for the over-all reaction. 

The concentration of permasurate ion is calcu- 
lated as about 10-” mole/liter. The value as- 
sumed for the activity of the permasurate ion is 
not too critical, for the hydrogen ion enters the 
Nernst equation to the eighth power and con- 
tributed heavily to the logarithmic term. 


RT 
E=E°—— 
nF 


92 
log [MaO.-] 


0.05928 


—0.146 = £°+0.101+0.160, 
E°(Ma—Ma04-) = — 0.41 volt. 


[H*], 


Whereas this value is an approximation, it 
should be pointed out that even a thousand-fold 
error in the assignment of the permasurate 
activity produces a change of less than 0.03 volt 
in the final value. 

This value for the standard electrode potential 
compares favorably with the known values for 
the electrode processes involving perrhenate and 
permanganate. 


E°=—0.15 volt for the reaction: 
Re+4H,0= ReO,-+ 8Ht+7e. 


The E® for the comparable reaction involving 
manganese is not given in the literature but can 
easily be calculated from the known standard 
free energies of the ions.® It is found equal to 
—0.78 volt for the reaction: 

Thus, the standard electrode potential of Element 
43 falls between the standard electrode potentials 
of its border elements in group VII (b) and ap- 
proaches the potential of rhenium more closely 
than it does manganese. The standard free energy 
of the permasurate ion is calculated to be 
—1.6X10*® calories from the above data. The 
standard free energies of the perrhenate and 


*W. M. Latimer, Oxidation Potentials (Prentice-Hall, 
Inc., New York, 1938). 
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permanganate ions are — 2.028 X 105 calories and 
— 1.006 X 10° calories, respectively.® 


OTHER OBSERVATIONS 


In an effort to secure further information con- 
cerning the oxidation-reduction of Element 43, 
some studies of the action of reducing agents on 
the oxidized 43 were made. 

The experiments were simple. To a given 
volume of a solution of 43 in sulfuric acid was 
added the reducing agent whose effect was to be 
studied. The reducing agents used were: metallic 
zinc, tin, nickel, lead, copper, and mercury; 
concentrated hydrochloric acid ; and a 20-percent 
solution of stannous chloride in 1:1 hydrochloric 
acid. In the case of the metallic reducing agents 
the solutions were tested for residual activity. 
Where the other reducing agents were used, the 
solutions were tested chromatographically on 
columns of alumina to determine whether the 43 
was anionic or cationic. 

Chromatographic studies were also conducted 
on the solution of active 43 used in the experi- 
ments on electrodeposition and on active solu- 
tions prepared by dissolving electrolytic deposits 
in 0.1N hydrochloric, sulfuric, and nitric acids. 

The chromatographic technique of determining 
the nature of ions has as its basis the well known 
facts (1) that an acidified column of alumina 
adsorbs only anions and (2) that an untreated or 
basic alumina column adsorbs only cations. 

Short columns, approximately one inch in 
length, were prepared in glass droppers using 
Fischer Adsorption Alumina, 80-200 mesh. The 
“acid” columns were prepared according to Zech- 
meister ;!° the ‘‘basic’’ columns were untreated. 
These columns were employed with the following 
technique : 

A five-milliliter portion of the solution under 
study was added by pipette and allowed to filter 
through the alumina into a ten milliliter volu- 
metric flask. The column was then washed with 
water until the filtrate was up to volume. This 
was followed by the collection of a second ten- 
milliliter portion of washings. A count was taken 
on the original solution and on each portion of 
filtrate to determine whether or not adsorption 
had taken place. 

LL. Zechmeister and L. Cholnoky, Principles and 


Practice of Chromatography (John Wiley and Sons. Inc., 
New York, 1941), p. 314. 
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The results of the reduction studies were as 
follows: 


Metallic reducing agents. 


All of the metals studied caused a more or less 
quantitative decrease in the activity of the 
original solution. This must be attributed to one 
or more of three possibilities: (1) reduction to the 
metal, (2) reduction to an oxide, or (3) physical 
adsorption on the “reducing” agent. 


B. Chromatographic results. 


1. The solution of 43 that was used for the 
electrolysis studies showed complete adsorption 
on an acid column and no adsorption on a basic 
column. This establishes its anionic character. 

2. Solutions prepared by dissolving electro- 
deposits of 43 in very dilute hydrochloric, sulfuric, 
or nitric acid gave no evidence of adsorption on 
acid columns, but after neutralization were 
adsorbed on a basic column. The activity must be 
considered cationic. 

3. Stannous chloride changed the activity from 
an anion to a cation. Such a change can only be 
attributed to reduction. Electrodeposition studies 
on the reduced solution might yield interesting 
results. 

4. Hydrochloric acid produced a definite, al- 


PERCENT DEPOSITED VS. CATHODE POTENTIAL 
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though not quantitative, change in the chroma- 
tographic behavior of the activity. In all proba- 
bility it is a reducing agent for 43, as it is for 
rhenium and manganese. 

The authors wish to extend their thanks to 
Dr. C. W. Balke, of the Fansteel Metallurgical 
Company, North Chicago, Illinois, for the sam- 
ples of high purity molybdenum. 


CONCLUSIONS 


From the foregoing experiments, the following 
conclusions regarding the electrochemistry of 
Element 43 may be drawn: 

(1) The element may be obtained in an oxi- 
dized, anionic form treatment with nitric and 
sulfuric acids. 
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(2) The element may be deposited electro- 
lytically on a platinum cathode, with reduction 
to the metal very probably the principal process. 
No anodic deposition is observed. 

(3) The deposition potential for the metal 
from its anionic form is —0.146 volt vs. the 
normal hydrogen electrode in a solution of pH 
2.36. This corresponds to a standard electrode 
potential of —0.41 volt for the process: 


Ma+4H,0 = Ma0O,-+8Ht+7e. 


(4) The anionic form of the element is reduced 
and removed from solution by the metals: zinc, 
iron, nickel, tin, lead, copper, and mercury. 

(5) Hydrochloric acid and stannous chloride 
reduce the anion to a cation. 
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Distribution functions, Fa(z, {n}), for multicomponent 
systems are defined proportional to the probability density 
that n molecules in an infinite isothermal system of 
fugacity set z will occupy the configurational coordinates 
symbolized by {n}. All thermodynamic functions may be 
obtained as certain sums of integrals of these distribution 
functions. These sums are always convergent, but im- 
practically slow in convergence for numerical use without 
further transformation. In particular, the grand-partition 
function, exp [VP(z)/kT], may be expanded in a power 
series in the fugacities z with coefficients given by integrals 
of the distribution functions Fn(0, {n}) at the fugacity 
set 0. As has been previously demonstrated for one com- 
ponent sysiems, this is shown to be a special case of a 
more general relation permitting the calculation of the 
distribution functions (and therefore the thermodynamic 
functions) for one fugacity set from those at another set. 
The function —kT In Fn(z, {n}) is the potential of average 
force, Wa(z, {n}), of mn molecules at the fugacity set z; at 
zero fugacities, the potential of average force is simply the 
potential energy. When only short range forces are present, 
by employing the procedure previously used for the im- 


1. INTRODUCTION 


HE sstatistical derivation of. the thermo- 
dynamic functions of a perfect gas* from 
* Publication assisted by the Ernest Kempton Adams 


Fund for Physical Research of Columbia University. 
{ Submitted in partial fulfillment of the requirements 


perfect gas, the integrals of the distribution functions 
may be simplified to integrals of lower order which become 
the coefficients in such expressions as the power-series 
expansion of the pressure increment in terms of the fugacity 
increments. These series have the advantage of a much 
higher rate of convergence than those involving integrals of 
the distribution functions, themselves; however, they do not 
converge when the initial and final fugacity sets bridge a 
phase transition. The general equations are applied to the 
imperfect multicomponent gas, to isotopic gas mixtures, 
and to condensed systems with the variables converted to 
the usual activities and activity coefficients. A change of 
independent variable with a corresponding change in 
coefficients then permits the evaluation of the pressure 
increment in terms of the solute concentrations. The 
equations for osmotic pressure are developed and found 
to be entirely analogous to those for the pressure of an 
imperfect gas. Finally, the osmotic pressure in a binary 
system and its behavior near phase transitions are dis- 
cussed with emphasis on the critical mixing phenomenon, 
for which the existence of a region of anomalous first-order 
transition is shown to be probable. 


a knowledge of the mechanical properties of the 
molecules has long been known. A simple state- 


for the degree of Doctor of Philosophy at Columbia 
University. 

*Here by “perfect gas’ we mean one which is non- 
gaan as well as having negligible intermolecular 
orces. 
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ment of the solution is as follows: Consider a 
single molecule of species s in a volume V and 
let €1, be its energy when in a particular quantum 
state Js. The partition function, Q,, for such a 
system is a function of T and V, and is defined as 


Q.(V, T) =D exp (—ex/kT). (1) 


Since the Hamiltonian is separable in the coor- 
dinates and momenta of the center of mass, the 
partition function (1) becomes a product of two 
factors. One factor, due to the motion of the 
center of mass, is, for macroscopic volumes, 
(2am,kT /h?)iV; the other, Qis, is due only to the 
internal quantum states: 


Q.(V, T) = is, (1’) 


where m, is the mass of the molecule of species s, 
and Qj, is the internal partition function, 


in which the summation runs over all internal 
quantum states j, of energy €js. 

The chemical potential per molecule, y,, of a 
perfect gas composed of NV, molecules of species 
sin volume V at temperature T is given by 


us =kT(In N,—In Q,). (2) 


If the particle density, p,=N./V, is large, a 
system composed of real molecules is no longer 
a perfect gas, and (2) ceases to apply. Neverthe- 
less, the formal method for obtaining the thermo- 
dynamic functions is still known: We now con- 
sider a system composed of the set N=MN,, 
No, «++, Ns, «++, Ne of molecules of species 
1, 2, ---, s, in a volume V. Let / bea 
running index representing single quantum 
states of this system, and E;, the energy of the 
Ith state. The partition function Q(N, V, 7) of 
such a system (which is calculated for the 
so-called canonical ensemble) is defined as! 


The work function, or Helmholtz free energy A, 
of the system is then given by 


A(N, V,T)=—-kTInQ(N,V,T). (4) 


1Cf. J. E. Mayer and M. Goeppert-Mayer, Statistical 
Mechanics (John Wiley and Sons, Inc., New York, 1940). 
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Two cases which have been successfully treated 
numerically by the use of this partition function 
(3) are those of a regular crystal (the Debye 
treatment which, however, is only approximate) 
and the perfect gas, which breaks down to Eqs. 
(1) and (2) since then, 


Q(N, V, T) (5) 


One may alternately employ the method of 
the grand-canonical ensemble, which will be 
elaborated upon in Section 4, since it forms the 
basis of the methods amplified here. This results 
in an expression for the pressure P in terms of 
the canonical gpartition function Q as defined 
above and the chemical potentials yu, of the 
various molecular species s: 


exp (PV/kT) 
[exp (Ls V,T), (6) 


where the summation runs over all sets of 
numbers N = No,---, ++, Ne of molecular 
species 1, 2, ---,s, +++, ¢ (O<N,< including 
N =0, for which we arbitrarily take 


QO, V, T) =1. 


Although Eqs. (3) and (4) do not appear to 
be fundamentally more complicated than (1) 
and (2), there is a quantitative difference which 
is so great as to become qualitative in character. 
It is usually possible to détermine the internal 
quantum states and the corresponding energies, 
€j, of the single-molecule system by a combina- 
tion of pure theory and spectrographic inves- 
tigation, even for a relatively complicated 
molecule. The summation involved in (1) may 
always be performed by a combination of 
analytical methods and straightforward, though 
frequently tedious, numerical evaluations. On 
the other hand, the quantum states and the 
corresponding energies, E;, of a system composed 
of many dependent molecules can only be evalu- 
ated in a few extremely simple cases, for (3) 
involves summation over some exp (10%) terms 
instead of the hundreds or possibly even only 
tens of terms entering in (1’’). Unless analytical 
simplifications are available, Eq. (3) is therefore 
purely formal and of no use as far as numerical 
approximation is concerned. 

Such analytical simplifications, based upon a 
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method first proposed by Ursell,? have been the 
subject of several recent papers.*-!° The dis- 
cussion of this article involves an attempt to 
present the mathematical procedures and argu- 
ments, particularly in connection with the 
method of the grand-partition function,!® in 
detailed form, and to generalize this method for 
multicomponent systems. Although naturally it 
is not to be expected that rigorous equations will 
ever be obtained which lead to ready numerical 
evaluation for the more complicated systems, it 
appears that certain relationships can be derived 
which have the rigorous generality of the Gibbs 


statistics. 


2a. SYMBOLISM 


In treating a problem involving the coor- 
dinates of an extremely large number of mole- 
cules it will be necessary to adopt a system of 
symbols which are sufficiently complex in nature 
to make the resulting equations appear relatively 
simple. Thus, while the meaning of a symbol may 
be at first difficult to bear in mind, the compact- 
ness of the final expressions more than com- 
pensates for this initial inconvenience. Accord- 
ingly, we make the following definitions: The 
symbol 

(7) 


indicates the configurational coordinates of the 
ith molecule of species s, and in general will 
include internal coordinates, zis, ***, aS 
well as the Cartesian coordinates, Xis, Vis, Zis, 
of the center of mass. 

The coordinates of m, molecules of a particular 
species s will be represented by 


{ns} =(1,), (2.), (is), (ns). (8) 


In general, any boldface symbol will indicate 
a set of numbers or variables, one for each species. 


(is) =Xisy Visy Zisy Qiisy ** * 


2H. D. Ursell, Proc. Camb. Phil. Soc. 23, 685 (1927). 

8 J. E. Mayer, J. Chem. Phys. 5, 67 (1937). 

4J. E. Mayer and P. G. Ackermann, J. Chem. Phys. 5, 
74 (1937). 


5 J. E. Mayer and S. F. Harrison, J. Chem. Phys. 6, 87. 


(1938). 

as 38) Born and K. Fuchs, Proc. Roy. Soc. A166, 391 
7B. Kahn and G. E. Uhlenbeck, Physica 5, 399 (1938). 

hem. Phys. 7, 1025 


2 (1941). 
10 J, E. Mayer, J. Chem. Phys. 10, 629 (1942). 
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Thus, as in Section 1, the bold-face n will be 
used to indicate a set of numbers of molecules 
of (a finite number, o, of) different species, n, 
of species s: | 


“Ne, (9) 


the number of molecules under consideration 
being* 


No, ***, Ns, ** 


n=> 


s=1 


(9’) 


The coordinates of all the molecules of such a 
set n will be represented by 


{mn} {me}, {me}, (10) 


The differential volume elements are similarly 
defined by 


d{n,} =d(1,)d(2.)- -d(i,) 
and 


(11) 
(12) 


d{n} (13) 


and are understood to be calculated in the 
Cartesian space of the atoms composing the 
molecules, i.e., JLx/q] is the Jacobian of the 
transformation from the generalized internal 
coordinates, Qiis, Geis, 

The ensemble or system treated will always 
be understood to have a fixed temperature 7. 
Its thermodynamic state might be fully deter- 
mined by giving T and the density set . 


O=P1, P2, ***, Ps, Pay (14) 


where p, is the number of molecules of species s 
per unit volume, 


Ps >= V, (15) 
or by specifying T and the chemical potential set 
U=Hi1, M2, ***, Msy Moy (16) 


with yu, the chemical potential of molecules of 
species s. Instead we shall find it more con- 
venient to use J and the fugacity set 


Z=21, 22, Bay (17) 


in which z,, the fugacity of molecules of type s, 


” Unless otherwise stated, sums and products over the 
species s will be assumed to have the limits 1<¢ s<co. 


| 
| 
(1939). 
-_ ® J. E. Mayer and Elliott Montroll, J. Chem. Phys. 9, ————— 
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is defined as proportional to exp (u./kT) with 
the proportionality constant so chosen that at 
infinite dilution z,=p,. This definition may be 
written 


,=[exp exp (—Hos/RT) ], (18)* 
where Pos is the density in molecules of species s 
per unit volume at the chemical potential set wo. 
The proportionality constant in (18) may be 
evaluated directly by considering a system of 
fixed volume V but containing a variable num- 
ber, N., of molecules, i.e., having a variable 
density set oo. By Eq. (2) and (15), we have 


lim [(Noe/Qz) exp (—Hos/RT) ]=1 


Qo—0 
= ( V/Qs) lim [Pos exp (—pos/kT)], 
Qo—0 


since V and Q, are independent of ,. (18) then 
becomes 
exp (us/RT). (18’) 


The concentration activity coefficient y,(z) of 
species s may be defined by the equation 


=2s/ps(Z). (19) 


The product of this dimensionless quantity with 
the concentration (density) of that species gives 
the concentration of an ideal substance having 
the same free energy. 

Since the sum 


occurs frequently it is convenient to introduce 
the notation 


2b. THE DISTRIBUTION FUNCTIONS 


The discussion of this article concerns itself 
primarily with certain distribution functions, 
Fa(z, {n}), of the coordinates of a set n of 
molecules. These functions are defined by the 
statement that in an infinite system of fugacity 


set Z, 

CII. ps"*]Fn(z, {n})d{n} 
is the probability that a set n of molecules will 
be found at the coordinates {n} in the con- 


* The fugacity z, and its thermodynamic relationships 
are more fully discussed in Appendix A. 


figuration element d{n}. For such values of the 
coordinates {n} that all intermolecular distances 
are large, this probability, integrated over the 
internal coordinates of each molecule, must have 
the average value [J], p.”* per unit volume in 
the 3n-dimensional Cartesian space of the 
centers of mass of the molecules. The definition 
of the functions Fn is therefore such that their 
average value, under these conditions, is unity; 
i.€., 


[.(z) Fn(z, is the probability) 
density of the configuration {n} in an 
infinite system of fugacity set z, 


and -(21) 


The definition includes the case z=0 (that 
is z.=0 for all s) in which case, for a clas- 
sical system, Fn(0, {n}) is proportional to 
exp —[Un{n}/kT], with Un{n} the potential 
energy of the molecules as a function of their 
coordinates {n}. The general significance of 
—kT In Fa(z, {n}) as a potential of average 
force is discussed in Section 6. 

_It will be shown that a detailed knowledge of 
these distribution functions would make pos- 
sible the calculation of all the thermodynamic 
properties of the system. 

As we proceed it will be evidently expedient 
to introduce new variables and functions, but to 
preserve the continuity of the development we 
postpone these definitions until they become 
necessary. For reference a summary of symbols 
is appended. 


3. SUMMARY OF THE METHOD 


Having defined the principal functions and 
quantities employed, the method and results can 
be summarized briefly before delving into the 
detailed derivations. 

It is a simple consequence of the equations of 
the grand-partition function (Section 5) that 
exp VP(z)/kT may be written as a sum over 
all sets n of terms of the form . 


if J Fa(0, {n})d{n}. 
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Similarly it will be shown that one can weite 
(TI {N}) exp VP(z)/kT 


as a sum of terms over all values of n, the term 
for n being 


/m nif Fa+w(0, {n-+N})d{n} 


(Eq. (33)). We thus have an infinite set of 
equations, one for each N, for; the unknown 
functions F(z, {N}), each equation involving 
on the right-hand side an infinite series. There is 
no reason to expect that from such a set one can 
obtain explicit solutions for the integrands in 
these coefficients in terms of integrals of the func- 
tions Fn(z, {N}), but just that is the case! By an 
algebraic substitution, it is then possible to write 
for any fugacity set z+yy a series for each function 
Fn(z+y, {N}) permitting the calculation of 
these functions from the Fn+m(z, {N+M});iie., 
one may calculate the distribution functions at 
z+y from a knowledge of those at any other 
fugacity set z. 

Now the particular feature of these equations 
which is of interest is the following: In going 
from any fugacity set z to another, z+y, the 
series involve only the difference in fugacities, 
Voy Vey Yo, raised to appropriate 
powers, the coefficients being functions of the 
original fugacities. The particular case for which 
21) is not unique; in 
this case the difference-fugacity y, is the com- 
plete fugacity. 

Thus the whole method of the grand-partition 
function is generalized. The thermodynamic 
properties and the distribution functions at any 
fugacity are, in principle, calculable from those 
of any other fugacity, the case: density, pressure, 
and fugacity, equal zero being only one special 
case for which, incidentally, the logarithms of 
the distribution functions are (in a classical 
system) proportional to the true potential ener- 
gies. The equations for calculating the properties 
from another set of fugacities are essentially the 
same as for the special case usually treated, the 
difference of the fugacities alone entering. 

Although these series converge in the mathe- 
matical sense for all values of y,, they are not 
suitable for numerical evaluations, for the con- 
vergence is far too slow. However, it is possible 
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to transform them into other series which, under 
certain conditions, lend themselves more readily 
to numerical approximation. One makes use of 
the fact that the functions Fn(z, {N}), according 
to (21), are unity (in the Cartesian space of the 
centers of mass after integration over the in- 
ternal coordinates) for finite N and large V over 
the considerable portion of space where all the 
molecules are distant from each other. We write 
the functions Fn as unity plus a series of cor- 
rection terms which differ from zero only where 
pairs, triples, or in general certain sets of mole- 
cules are close together. By appropriate defini- 
tions of the correction terms the equations still 
involve no approximations. By this means the 
integral over {M} of each Fn+m may be written 
as a sum of products of integrals over these cor- 
rection terms involving fewer molecules than M. 
The coefficients occurring are such that the 
complete summation over M may be written as 
an exponential of a power series in the fugacity 
differences y and the activity coefficients y(z). 
From the first equation so obtained, 


[P(z+y) —P(z) ]/kT 
is expressed as a power series in y,/ys(Z) while 
the higher equations yield similar series (in the 
same variables) for the quantities 


In Fn(z+y, {N})—In Fr(z, {N}), 


the coefficients in the latter case being, of 
course, functions of the coordinates {N}. These 
power series have the property that for small 
values of the arguments they converge numeri- 
cally, in the sense that the terms decrease rapidly 
in magnitude, and may therefore be used for 
numerical evaluation. They are not, however, 
convergent for all values of y, and indeed diverge 
for such values of the temperature and fugacities 
that any phase transition occurs in the system 
between the fugacity z and z+y. They thus 
lead to equations which are valid for the calcu- 
lation of the properties at the fugacity set z+y 
from those at z even if the phase present at z is 
condensed, provided only that there is no phase 
change between the two fugacities z and z+. 
The application of the general equations to 
binary liquid mixtures leads to the laws of ideal 
solutions and the specifications which must be 
met in order that the ideality be preserved 
throughout the whole concentration range. 
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4. THE GRAND-PARTITION FUNCTION FOR MULTICOMPONENT SYSTEMS 


For the sake of completeness and consistency of notation we devote this section to a short deriva- 
tion of the equations for the grand-partition function ' '? for an ensemble of constant energy having 
maximum-entropy. 

We consider a grand canonical ensemble consisting of a total volume MV divided into M similar 
regions each of volume V, the total energy and the total number of molecules of each species of the 
ensemble being fixed. The regions are separated, in thought, by semi-permeable membranes through 
which molecules and energy may flow, the separation being such that the total-phase space is the 
product of the independent phase spaces of the regions, i.e., the molecules in different regions 
interact with each other only negligibly. 

Consider a number, w, of quantum states having energy between E and E+AE, available to a 
system composed of the set N of molecules in one such region.* Further, let m be the number of such 
systems in these w states and let Q(m) be the total number of quantum states available to the 
ensemble of systems with the inhibition that exactly m are contained in the w states. Since the total 
number of particles in the ensemble is constant, as is the total energy, the criterion of equilibrium 
is that the entropy, S, of the ensemble be a maximum. Now the number of distinct ways in which 
we can select m systems is M!/m!(M-—m)!, and one thus obtains for 2(m) the expression: 


Q(m) (22) 


where 2,(4—m) is the total number of quantum states available to the remaining M—m systems. 
Employing a familiar definition of the entropy, 


S(m) =k In Q(m) =k In M!—& In m!—k In (M—m)!+km In o+ S,(M—m) 
and 
S(m—1)=k In M!—k In (m—1)!—k In (M—m-+1)!+(m—1) In o+ S,(M—m-+1). 


At maximum, 


AS(m) =S(m) —S(m—1) =0 


(M—m+1)o 
=k In +AS,(M—m) =0, (23) 


where AS,(M—m) =S,(M—m)—S,(M—m-+1). 

Since the change from m—1 to m of the number of regions in the w quantum states involves in the 
rest of the ensemble a decrease of V in the volume, a decrease of NV, in the number of molecules of 
type s, and a decrease of E in the energy, 


OS, OS, OS, 
AS,(M—m) = v( ) -> ) ) 
OV/JNE ON,/ Ni, V,E OE J 


(24) 


1 R. C. Tolman, Phys. Rev. 57, 1160 (1940). 

®E. A. Guggenheim, J. Chem. Phys. 7, 103 (1939). 

* In the total-phase space the component parts belonging to the different regions are numbered due to the fact that 
these regions occupy different parts of the physical three-dimensional space, so that although for every quantum state 
of one region corresponding to a fixed number of molecules in that’ region in a definite state, there corresponds a similar 
quantum state of every other region differing only in a spatial translation, the corresponding quantum states of different 
regions are not identical. 
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in which it is understood that M approaches infinity. Thus, as M@>>m>1, (23) and (24) give 


m 
(—PV+N-u—E,))/kT (25) 


@ 


=@(INV, 


the probability that a randomly selected region of volume V contain exactly the set N of molecules 
and occupy a particular one of the w quantum states of energy E;. The sum of @unv over all possible 
quantum states / of each N, and over all numbers of molecules N, is unity: 


> =1=exp (—PV/kT) exp (N-u/kT) exp (—E,/kT), (26) 
N 1 N>O l 


where £, is the energy of the /th quantum state. We thus obtain for the grand partition function of 
a multicomponent system: 


exp (PV/kT) = exp (N-y/kT) exp (—E,/kT) 


N>o 


(27) 
= [exp (N-y/kT) JQ(N, V, T), 
N>0 


in which Q(N, V, 7), Eq. (3), is the partition function for the canonical ensemble containing the 
set N of molecules in the volume V at temperature 7, and has the value unity for N =0. 


5. GENERALIZATION OF THE GRAND-PARTITION FUNCTION 


In the last section it was shown that the probability that there be exactly the set N of molecules 
in the volume V and that this system be in the quantum state 7 is: 


@mnv =exp (25) 


where £; is the energy corresponding to the quantum state /. The probability that these molecules, 
being now for the moment numbered and collectively in the quantum state /, will occupy (re- 
spectively) the coordinates of configuration {N} within the volume element d{N} is |y.{N}|?d{N}, 
where y; is the normalized wave function of the system. 

Let @nv{N} be defined by the statement that Pnv{N}d{N} is the probability that there is exactly 
the set N of molecules in V and that these are located at the coordinates of configuration {N} 
within d{N} irrespective of the quantum state of the system. From (25) one sees that 


Onv{N} = {exp (1/kT)[—PV+N-y]} exp (—Ei/kT)| N} (28)* 


If the volume V is known to contain the set N+n of molecules then the probability density that 
there be N molecules at the positions {N} irrespective of the positions of the other n molecules in 


the volume V is 
SSS v{N-+nj}d{n}, (29) 


since any N, of the N,+m, molecules of species s may fill the {N;} positions, and having chosen a 
particular set NV, of molecules of each type s, all permutations of the molecules among the positions 


* This is the probability that, having numbered the molecules, the ith molecule will be at the ith position, etc., and 
thus if integrated over all values of the coordinates within V, gives the probability that exactly this number of molecules 
_ are in V regardless of position. The probability that the positions {N} be occupied by molecules of the appropriate 
—_ ae of which molecule occupies which coordinate set, is (I1,V,!) times (28); this latter quantity was used 
in reference 10. 
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are allowable. If we now sum (29) over all values of n, we get the probability density that the con- 
figuration {N} is occupied by the appropriate species of molecules at this temperature T and fugacity 
set Z, namely 


(N.+n,)! 


Ng: 


(1 (—PV/AT) exp ((N+0)-w/#7) I] 


XE exp (— f f J (30)* 


in which (28) has been substituted for @N+n,v, and the expression defined in (21) has been used. 
Employing our previous definition of the fugacity (18’) in (30) we obtain 


|| \?exp (—E,/kT)}d{n}. (31) 


Taking the limit of this equation as all z, approach zero, the only finite term remaining in the sum 
on the right is that for which n=0; therefore, 


V 
“ 8 
and (31) may be written: 


{N}) =exp (—PV/kT)X Frist, {N+n})d{n}, (33) 


Zs n,! 


or, in more symmetrical form, introducing P(0) =0, y,(0) =1, 


exp (VP(2)/AT)| TIC.(2) Ful, IN})= | f f f exp (VP(0)/kT) 
nz V 


n,! 
Employing (21) for N=0, Fo(z, {0}) =1, and we have, as a special case of (34), 


2," 
exp (VP(2)/ET) =¥ | TI fe {n})d{n}. (35) 
Now form the sum 


m20 m, | 


* Actually, an equation of this type defines FN(0, {N}) only for such values of the coordinates {N} of the molecules 
as lie within the volume V for which the quantization of the right-hand side is performed. However, such a function 
defined within V is asymptotically equal, at al! values of the coordinates {N} for which it is defined, to the function 
defined by (30) for any volume greater than V which includes V provided only that V is large enough and of such a 
— that surface effects may be neglected. We may consequently regard the distribution functions as defined over 
infinite s . 

Slater, Phys. Rev. 38, 237 (1931). 

In reference 10 the quantities {[y(z)]-"} F,(z, were abbreviated as Q,(z, {m}). 
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Substituting (34) for the integrand, (36) becomes 
)magmetns 
-———) f f- J Frimin(0, {N-+m-+n})d{m-+n}. 


Regrouping the terms according to the value of m+n=j, we have 


=Fn(0, {N}). 


Therefore, (36) may be written 


X Fn+m(z, {N+m})d{m}. 


(37) 


By rewriting (34) for the fugacity set z+y, and employing (37) in the integrand we obtain 


anol 
exp TI {N}) = Ul ( Ys) ) 


n2>0 m2>0 Ms '‘n,! 


Regrouping terms according to the sum n+m=M, and again using the multinomial expansion the- 
orem, we find 


xexp J | (N+M})d{M}. (38) 


This is a generalized form of the grand-partition function for a multicomponent system, enabling 
the calculation of the pressure, P(z+y), the activity coefficients y,(z+y), and the distribution func- 
tions, Fn(z+y, {N}), at any fugacity set z+y, in terms of the pressure, P(z), the activity coef- 
ficients, y,(z), and the distribution functions, Fn(z, {N}), at any other value z of the fugacity set. 
Thus Eqs. (34) and (35) are special cases of (38) with z=0, while (37) is the special case y, = —2s, 
y+z=o0. 

So far the equations derived have not been restricted to any particular kind of system. As will be 
shown immediately, all of the series converge because of the finite size and low compressibility of 
real molecules, this being equally true even if, for instance, the fugacity sets z and z+-y in Eq. (38) 
correspond to two different phases.* However, as has been previously pointed out (cf. Section 3), 
the convergence of these series is too slow to permit, even in principle, numerical evaluation. The 
demonstration that the series are absolutely convergent is simple: the functions Fn+m(z, {N-+M}) 


* Note that, owing to the omission of gravitational forces, a condensed phase will not necessarily be found in any 
particular region of the containing vessel. 


7) 
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are all positive, and for exceedingly large values of M approach zero in value for all values of the 


coordinates {N-+M} within a finite volume V, for all real systems for which the molecules repel 
each other at small distances. Since, then, the integrals approach zero in value for infinite M values, 
the convergence is assured. However, the term which gives the largest contribution to the sum is 
that for which N+ M represents the average molecular population of the volume V, and for ordinary 
macroscopic volumes is therefore of the order of magnitude of 107°. 

An enormous improvement in the rate of convergence for systems satisfying certain requirements 
may be obtained by employing a method first suggested by Ursell,? and extended by Kahn and ~ 
Uhlenbeck’ to quantum mechanical systems. This procedure will be discussed in Section 7, but in 
order that the development there can proceed uninterruptedly, we pause to introduce the concept 
of the potential of average force which will be needed in simplifying later equations. 


6. THE POTENTIALS OF AVERAGE FORCE 


If we may treat the system consisting of N molecules in volume V semi-classically, i.e., if the 
spacing of the quantum mechanical energy levels is small in comparison with kT, the probability 
density of finding the system in the configuration {N} is given by the expression 


Semi classical 


1 


1 
WQ TL. Ne! 


where f is the total number of degrees of freedom of the system, and the classical Hamiltonian, 
H({N}, {Pw}), is the sum of the kinetic and potential energies which in turn are functions of only 
the momenta and coordinates, respectively : 


H({N}, {Pw})=7{Pn} + Un{N}. (40) 


f fee (Pw})/kT]d{Pw} (39)* 


By performing the integration indicated in (39) and recalling the equation obtained in the previous 
section for the distribution functions at zero fugacity, 


Fx(0, {N}) -[n(5) ™ exp (32) 


it is evident that there is a direct proportionality between Fn(0, {N}) and exp [— Un{N}/kT] for 
any given N. Indeed, by choosing the zero of potential energy correctly one may write: 


Fn(0, {N}) =exp [— Un{N}/kT ]. (41) 
The correct zero of the potential energy may be easily obtained, but it is unessential in the argument 
which follows. (41) must, of course, satisfy the normalization condition (21). 


If, for convenience of notation, we now consider a system of n+N molecules in volume V, the 
force acting along the coordinate gq of the molecule 7, occupying one of the {N,} positions is given by 


(42) 


fais{n +N} =—- 


A comparison of Eqs. (28) and (32) shows that the probability density that there be exactly the 
set n+N of molecules in V and that of these, n be located at the coordinates of configuration {n}, 


* See, for example, reference 1, p. 231. 
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‘and the other N be distributed in any way among the {N} remaining coordinates is (see Eq. (29)) 
(n.+N,)! 


Ns! 


nst+Ne 


poss vin+N} =exp (— TI {n+N}). (43) 


n,! 


The average force, (feis{N}), averaged over all positions of the coordinates {n}, may now be found by 
multiplying (43) by (42), integrating over the coordinates {n}, summing over all values of n, and 
dividing this result by the sum and integrals of (43); that is 


exp (—VP(2)/kT) f f f (-—*) exp (—Ua+n{n-+N} /kT)d{n} 


(fais{ N}) = ; 
or, 
0 


in which (33) has been used. The quantity —kT In Fn(z, {N}) is therefore such that its derivative 
with respect to one of the coordinates {N} gives the statistical average force, as a function of {N}, 
along that coordinate on the molecule occupying that position, the weighted average being taken 
over all positions {n} of the other molecules in the system; or what is equivalent, the average is 
made over a long period of time during which the {N}-coordinates are fixed and occupied by an 
appropriate set of molecules, and the remaining n-molecules in the system are allowed to move 
freely. It thus is the potential of average force, denoted by Wn(z, {N}), and we write 


Wn(z, {N}) =—T In Fn(z, {N}). (45) 
In the limit of zero fugacity, it is evident from Eq. (41) that 
Wn(0, {N}) =Un{N}. (46) 


In previous treatments it has been convenient to consider the potential energy Un{N} as a sum 
of component potential energies, ~;;(r;;), due to the interaction between pairs of molecules only. 
An extension of this method to the potentials of average force leads to some simplification of the 
equations of Section 7. If we denote a subset (consisting of the coordinates of n molecules) of the set 
of coordinates {N} of N-molecules by the symbol {n}n, we may express the function Wn(z, {N}) in 
terms of the component potentials of average force by the equation 


Ww(z, {N})=2 {njnwa(z, (47) 


in which the sum is taken over all possible subsets, and the terms wa for n>2 may be thought of as 
corrections for the deviation of Wn from a sum of pair terms only. The inverse of (47) is given by 


wa(z, {n}) = 27 (48) 


which is proved in the Appendix. It has been general practice to neglect terms wn for > 2. 

An interesting and important observation regarding the general Eq. (33) for the distribution 
functions is that the integration over the internal coordinates of the N molecules may be performed 
before integration over the center of mass coordinates, {N,}, giving a distribution function de- 
pendent upon positions but independent of the internal coordinates of the molecules. This new 
distribution function, which will be referred to hereafter as the “spatial distribution function,” 
may then be treated just as though it were for a system of monatomic molecules, except that the 
quantities corresponding to the potentials of average force are no longer even approximately equal 
to the sums of pair terms. 


29)) 
43) 
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7. EXPANSION AND INTEGRATION OF THE DISTRIBUTION FUNCTIONS 


According to (21), the spatial-distribution functions, Fm(z, {M.}), are unity for finite M and large 
V in the considerable portion of configuration space where all intermolecular distances between the 
M molecules, are large. If it happens that the force between groups of molecules drops rapidly to 
zero at separation distances comparable to the molecular size, we may write a given spatial-distribu- 
tion function as unity plus a series of correction terms which differ from zero only where pairs, 
triples, or in general certain sets of molecules are near each other. It will be shown convenient to 
choose the form of these correction terms so that Fm(z, {M,}) may be written as a sum of products 
of functions each of a smaller number of coordinates (Eq. (49)). The procedure is formal and 
generally rigorous, but to facilitate visualizing the method we apply it here to the spatial-distribution 
function only, and at such a fugacity set z that the system is gaseous or liquid; actually we later 
deal (Eqs. (49) and (50)) with the complete distribution function (i.e., including internal coor- 
dinates) for whose expansion an analogous explanation may be given. 

Before proceeding to an explicit definition of the correction terms it may be well to stress this 
feature of the distribution functions: When the center of mass coordinates, {M,} admit of a natural 
spatial separation into several mutually exclusive subsets of neighboring molecules—that is, when 
the M molecules may be segregated according to position into several groups such that all the 
molecules of each group are at large distances (compared with the molecular diameter) from all 
the molecules of every other group, yet are ‘‘close’”’ to each other—the spatial-distribution function 
F(z, {M.}), being a probability, must be expressible as a product of probabilities, (i.e., distribution 
functions), one for each of the simultaneous independent events of having a particular subset of 
molecules close together. Thus, when one pair is in conjunction, with all other molecules mutually 
distant, we correct the value unity by adding a term for this pair. However, in order that the above 
requirement be fulfilled when a second pair of molecules are near each other we must add, in addition 
to the correction term for the second pair, another term which accounts for there being simul- 
taneously two pairs present. This additional term is readily seen to be the product of the two different 
pair-correction terms. 

If a triple were present we would have, in addition to the term corresponding to the triple, terms 
corresponding to all possible pairs which may be formed from the triple. In other words, if the coor- 
dinates {M,} are such that a particular aggregation exists, the expansion of the distribution function 
involves not only the term corresponding to this particular aggregation but in addition, correction 
terms for all lower aggregations which can be obtained from this aggregation by a process of dis- 
sociation. We begin with all molecules far apart, and correct first for pairs, then for triples, etc., 
until the desired state of aggregation is attained. 

The correction term for having simultaneously pairs, triples, etc., must vanish if any one of the 
groups for which it corrects is dissociated. A simple form for the correction term, and one which 
satisfies the multiplicative property of the distribution functions discussed above, is thus a product 
of factors, each factor corresponding to a particular group and having the property of vanishing 
when that group dissociates in any way. 

We denote as before (the coordinates of) a subset of m particular molecules of the set M by the 
symbol {m}m, and a complete set of k unconnected (i.e., mutually exclusive) subsets {m;}m(1<i<k) 
of the set M by the symbol {k{m;}m},, with 


k 
> m;=M. 
i=1 
We note that two subsets of identical size are different if they do not contain the same particular 
molecules. Let the correction term at fugacity z for {m;}m proximate molecules be gm(z, {m,}™m) ; 
then the correction term for having the set of subsets {k{m;}m}, will be 
k 


II gm(z, {m,}m), 


i=1 


287 
44) 
‘ive 
N}, 
ken 
2 is 
an 
ove 
45) 
46) 
um 
ly. 
the 
set 
) in 
47) 
as 
by 
48) 
ion 
ed 
le- 
ew 
he 
lal 


288 W. G. McMILLAN AND J. E. MAYER 


and we may write the expansion of Fm(z, {M}), defining the g-functions implicitly, as 
k 
Fu(z, IT gm(z, {m,;}™), (49) 


where the sum is taken over all possible sets of unconnected subsets. An explicit definition of the 
g-functions is given by the inverse of (49): 


1 
gm(z, {m}) => {/{M;}m} II Fu(z, {M;}m), (50) 
which is proved in the Appendix. Note that, according to (49), if all the positions {M} are widely 
separated the integral of Fm over the internal coordinates has the average value unity, so that the 
average value of g,(z, (z)) is unity when integrated over the internal coordinates of the ith molecule. 
If the fugacity set z has such a value that the system is liquid or gaseous then the integral of g:(z, (z)) 
over the internal coordinates is unity. If the system is crystalline the integral over the internal coor- 
dinates may lead to a triply periodic function of the coordinates of the center of mass,!° but the 
average value over an appreciable range of these coordinates is still unity. 

The functions gm(z, {m}) are so defined that they approach zero if the values of the coordinates 
{m} are such that there exist two subsets {m’}m and {m}—{m/’}m with all distances large between 
any molecule of {m’}m and any molecule,of {m}— {m’}m. The proof that this is the case is now given. 
By the multiplicative property of the distribution functions the following statement holds for Fu: 
if the coordinates {M} are such that there are two mutually exclusive subsets {m’}mand {M}—{m’}m 
with all distances between the two subsets large, then 


F(z, {M}) =Fm-m’(z, {M—m’}m) Fm’(z, {m’}m). 


This relation is evidently consistent with (49) for Fm, Fm-—m’ and Fm’ only if all functions gm of the 
coordinates of sets {m}m containing members from both {m’}m and {M—m’}m are zero. 

Since (38) necessitates the integration of the function /n+m(z, {N+M}) over the coordinates 
{M} we require an expansion of this function similar to (49). If the coordinates {N-++M} are such 
that they fall naturally into two isolated groups of N and M molecules respectively at the coordinate 
sets {N} and {M}, we may write 


Fu+m(z, {N+M}) =Fwn(z, {N})Fm(z, {M}) 
k 
=Fr(z, {N}) tk{mi}m} II gm(z, {m;}™M). 
However, if it happens that some of the N molecules are close to some of the m we need an additional 
correction. Denoting a subset of (the coordinates of) n particular molecules of the set N by the 


symbol {n}n, the necessary corrections, gnm(z, {n}n, {m}m), may be applied to the functions gm in 


the form 
gm(Z, {m;}) +> {n}N gam(z, {njn, {m,}™), 


and our complete expansion for Fn+m becomes* 


Fn+m(z, {N+M}) =Fn(z, {N}) {k{m;}m}, I [gm(z, {m;}m)+ {n}N gam(z, {n}n, {m;}M) ]. (51) 


If we now integrate this expression over all values of the coordinates {M}, the integrals 
= ff f {mn} (52) 
Vv 


*It should be noted that, though we wish3to neglect higher terms in the numerical evaluation of some of our later 
expressions, these complete expansions involve no approximations whatever. 
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(IL. m,')bam(z, {11} x) = f f f gam(z, {m1}, {m} {m} (53) 


depend only upon the numerical value of m =m, mez, +++, Ms, +++, Me, and not upon which of the 
M, molecules constitute the m,. (Note that if the interaction forces, i.e., the correction terms, drop 
sufficiently rapidly to zero with increasing separation distance of one molecule from all the other 
molecules of a group, the cluster integrals, bm(z) and bam(z), of the first and second kinds defined by 
(52) and (53) are independent of the volume of the system.) A subdivision of the {M} coordinates 
into um subsets of size m may be obtained in J], M.!/[]m([]. ms!)"um! different ways, all of which 
appear in the sum (41). Integration over {M} thus leads to the equations 


ff- {N+M})d{M} 


= Fu(z, (N}) if f- J Cem(z, {mn;} mt) +30 gam(z, Jd {m,} 


! 


in which the sum over um is subject to the restriction }°m mum=M, i.e., m m,um= M, for all s. 
Now employing (54) in (38) we obtain 


exp (VP(z+y)/kT)Fn(z+y, =exp (VP(z)/kT) Fn(z, 


M>0 “m m 


The summation over M removes the restriction on the ym, and it is seen that the sum in the right 
of (55) is the expansion of an exponential, so that 


exp (VP(z+y)/kT)Fu(z+y, {N})/II =exp (VP(z)/kT) F(z, {N})/TI. 


m20 


¥s(Z) 


for all values of the fugacity set y which are smaller than the corresponding set of associated con- 
vergence radii of the series in the last expression. If we now employ (45) with the expansion (47) 
for the potentials of average force, and collect terms in (56) according to their dependence on V 
and the coordinate sets {n}n, we find* 


P(z+y)—P(z)= eT TE (57) 


m20 


Wi(Z+Y, (t.)) +RT In =W1.(Z, In dis, m(z, I [y./ve(z) (58) 


m20 


and 
wna(z+y, {n}w) =wn(z, {n}n) —kT 2 bam(z, {n}w) II (59) 
mz 8 
where the subscript 1,,m of 6 in the sum of (58) means that all m, are zero except n,=1. 
These equations are analogous to those obtained by Mayer,’ and reduce to those for a one 
component system. Particularly, equation (57) permits the calculation of the pressure at the fugacity 


* For convenience in writing summation limits, bo is taken equal to zero. 
{ Incidentally, Eqs. (43) and (44) of this reference contain several incorrect signs. 
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set Z+y from the cluster integrals and the activity coefficients evaluated at the fugacity set z, 
whenever the sum converges. This remark will be elaborated upon in Sections 8 and 9. 


8. THE CHANGE OF INDEPENDENT VARIABLE 


- In the integration of the g-functions (Eq. (52)) for simplified systems in which wa=0 for n>2, 
it has been found convenient to use irreducible integrals, Bn, which may be physically defined for 
such systems as the integral of a certain function of the potential energies between pairs of molecules 
in a group which is more than singly connected, i.e., one in which it is possible to pass from any 
molecule to any other molecule of the group by at least two alternate paths having only the initial 
and final molecules in common. The cluster integrals, bm(z), have been found to be sums of products 
of these irreducible integrals." 

The physical meaning of the irreducible integrals is somewhat obscured if the molecules are not 
perfectly symmetrical or if, in general, there are internal coordinates. However, as was previously 
pointed out in Section 6, an examination of Eq. (33) shows that we might have integrated Fn(z, {N}) 
immediately over the internal coordinates and applied the discussion of Section 7 entirely to the 
spatial distribution functions so obtained. Under these conditions the irreducible integrals have a 
perfectly definite physical meaning as solutions, in terms of the bm(z), of the implicit equation which 
defines them (see Eq. (60) below). In the actual calculation of the cluster integrals one would normally 
perform the integrations leading to the irreducible integrals first; add to this the simplification of 
having the readily measurable density rather than the fugacity as independent variable in our 
equations when the irreducible integrals are employed and there results sufficient reason for their use. 

The expression of the cluster integrals in terms of the irreducible integrals has been solved for one 
and two component systems by Mayer,®" and for any number of components by Fuchs." A brief 
outline of the general mathematical theorems proved by Fuchs will be given, and these will then be 
used to convert our series for the pressure (57) into one having irreducible integrals as coefficients. 

The solution to the combinatorial problem of expressing the cluster integrals, bm, in terms of the 
irreducible integrals, Bn, for a system of ¢ components may be written 

{m.n,Bn}’sn 
m,bm = (m,— 5st) vin}, Il II (60) 


os=] 0D vsn! 


in which the summation extends over all integers vn 20 subject to the condition that 


t=1 2 


and | |, symbolizes a determinant of oth order. Fuchs further shows that (60) may be written as a 
Cauchy integral: 


where G(é, Bn) and its derivatives are defined by the equations 
Bu) = Ba] 
n20 8 


dG 
G,(&, Ba) -(—) 
E. Phys. Chem. 43, 71 (1939). 


Fuchs, Proc. Roy. Soc. A179, 408 (1942). 
: For convenience Bo is taken equal to zero. 


and 
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Multiplying (61) by J]. Z."* and summing over m2 0, one obtains, on integration, the equations 
Z,GA(Z, bm)=Y,, r=1,2,°++,¢, (63) 


where the Y, are solutions of the equations 


Z,=Y,exp [[—G,(Y, Bn)], (64) 


Employing (63) and (64) it is then easy to show that 


G(Z, bm) =G(Y, Ba) +E Y[i— —G, (Y, Bn) }. (65) 


8=1 


This last result (65) is precisely what is needed to enable us to express the pressure difference 
between the fugacities z+y and z in terms of the irreducible integrals, for (57) may be written in 
the form 


m 20 Zs 


with Z,=¥.p.(Z)/Z.. It is evident that when 2, is held fixed, 


_& 
T, Zs Ys 


RT aZ, kT 


or 
+ y) 


(67) 


(66) with (67) in (65) then gives 


s=1 n20 s=1 


Since the appearance of a new phase in general causes singularities in the thermodynamic func- 
tions such as that given by (66), we conversely associate singularities in the fugacity function, 
G(Z, bm), with phase transitions. By the general theorem for implicit functions,’ the o-Eqs. (64) 
inversely define the variables Y in terms of the variables Z whenever the Jacobian of the trans- 
formation, J[Z/Y], does not vanish or become singular anywhere in the accessible domains of the 
variables. This Jacobian, by (64) has the value 


Z 
=| Y.G.r(Y, Bn)|, exp [x G.(Y, (69) 


s=1 


from which it is evident that a singularity in the function G(Z, bm) may arise either from singularities 
in the functions G(Y, Bn) or from the vanishing of the determinant in (69). This remark will be 
employed in some of the arguments of Section 10 in applying (68) to solutions. 


9. RESULTS 


The net results of importance for the evaluation of the pressure and, therefore, through the use 
of the equations of appendix A of all of the more usual thermodynamics functions, may now be 
summarized. 


va 1 for eam Goursat-Hedrick, Mathematical Analysis, first edition (Ginn and Company, New York, 1904), 
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Equation (57) expresses the pressure change in the system as a power series in the fugacity change: 


P(z+y)—P(z) =kT bm(z) (70) 


in which it is important to note that the coefficients of all the linear terms in the sum are unity, that 
is, bm(z) =1 for the subscript set m, =4,5. 
Using the thermodynamic relation (A13) of the appendix, 


LaP(z+y)/dy. Jr, Z, yr =kTp.(Z+y)/ (2+ Ys) 
applying this to (70) one obtains 


p(z+y) = mbm(2z) Il (71) 


Ys mz20 


as an expression for the density in molecules per unit volume of species s. 
Equation (64) with Z,=[y,/y,(z) ] and Y,=[y,/y-(z+y) ], Eq. (67), leads to the relationship: 


In y4(2) —In y.(2-+y) =~ Ly (72) 
while (68) is 
Vs 


P(a-+y)—P(2) =kT\ aa 
8 n20 8 ys(z+y) 
If the species of molecules were such that.all the mutual potential energies were independent of 
the species (e.g., for isotopes which are distinguishable in principle but negligibly different in me- 
chanical properties) then 


bm 


so that the integrals depend only upon the number of particles and not on their kind. If, in addition 
z=). 2,and p=). ps, the activity coefficients for all species are identical and the foregoing equations 
become 


where, as usual m=)°, m,, and 


m m! Vs ms 


m>0 m>o [],m,! 


() 


m20 8 m20 


In y(z) —In y(s+y) |", (72’) 


y 


where (71’) and (72’) are obtained from (71) and (72) by a summation over s. These equations are 
seen to be just those for the one component case.!° 


and 
| 
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The application to a gas may be made most conveniently by setting z=0, in which case y,(z) =1 
for all components, and y,/y.(z+y) while P(z) =0, so that Eqs. (70) 
to (73) become 

P(y) =kT bm II 


m2>0 8 


yet bm II 


m>1 8 


n>0 r 


p(y) = IT (75) 


1 
In Y,=In p.(y)-— > II (76) 
Ps n>0 


where the superscript zeros are used on the integrals to indicate that they are evaluated at zero 
fugacity of all components. The second line of Eq. (74) is written in this form to emphasize the fact 
that the coefficients of all the linear terms are unity. 

Essentially these equations for gaseous mixtures have been given before by Fuchs'® and their 
equivalents for two components were derived by Mayer.'* They deserve no further discussion here 
but are listed for comparison. 

The more novel application of Eqs. (70) to (73) is given if we consider the case of one or more 
solutes in a solvent.* For this case we shall find it convenient to use as the reference fugacity set z, 
at which the integrals are calculated, that at which the fugacities of the solutes are zero and the 
fugacity of the solvent is such that the pure liquid has some convenient pressure P*, say one 
atmosphere, or possibly its orthobaric vapor pressure. These conditions will be referred to hereafter 
as the standard state, z*. The activity of each component will be chosen with such a proportionality 
factor that it becomes equal to the density or, in more usual chemical terminology, the concentration, 
of that species in molecules per unit volume as the composition approaches that of the standard state. 

The notation will be altered to emphasize the unique position of the solvent, for which Greek 
symbols will be employed, while Latin letters will be adopted for the solutes. The solute quantities 
which appear turn out to be their concentrations in molecules per unit volume (which we shall now 
write ¢, instead of p,) and their activities a, in the same concentration units. The solvent quantities 
that enter are its concentration p in molecules per unit volume, its activity a (equal to the concen- 
tration in the limiting case of the pure solvent), and the activity coefficient y =a/p. For the standard 
state these quantities will be marked with an asterisk. 

If in our previous notation we choose the subscript zero to indicate the solvent} the standard state 
z* is given by 2,=2.6o;, and the new symbols are related to the old through the equations: 


p=p(z*+y) (77a) 
p* =p.(z*) =a* (b) 

a = (c) 
| Y0(z*+y)/y0(z*). (d) 


* J. G. Kirkwood has treated similar problems by a onan method. See J. Chem. Phys. 3, 300 (1935) and Chem. 
Rev. 19, 275 (1936). 
Tt So that the solutes are numbered from 1 to ¢. 


Solvent quantities + 
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p.(z*+y) (2,* +.) (78a) 
Solute quantities } 0=p,(z*) (b) 
ds =y,/7.(z*). (c) 


Since only the quantity y,/yo(z*) and not a@ itself appears in the equations we assign it the symbol 
Aa, 


Aa =yo/yo(2*) =a—a*, (79) 


it being the change in the solvent activity in’ going from the standard state z* to the state z*+y. 
For a solute, the change in activity between these two states is the total activity, the activity in 
the standard state being zero. The dissymmetry of the equations results entirely from the choice 
of the standard state. 

Employing the new symbols, Eqs. (70) to (73) becomeft 


P(a, a) =P*+kT du, m b%um(Aa)* ay"; (80) 


p(a, a) m a,""; (81a) 


Aa\’ 
Ina=In p—-— > II (b) 
Aa yn Y 
a) = m II- (82a) 
In a,=In > n,B*m II (b) 
Cs y,n 
Aa 
8 Y »n Y r 


where the asterisk on the integrals 6 and B. indicates that the integration is carried out with the 
potentials of average force occurring in the pure solvent, that is at infinite dilution of all solute 
species. 
Of particular interest is the case for which Aa =0, i.e., the activity of the solvent is held constant. 
The pressure difference 


Pla a) —~P(p*, 0) =n(p*, a), (84) 
is the osmotic pressure. We find 
m(p*, a) =kT b*om a,”", (85a) 
m>1 


where the second line is so written to draw attention to the fact that the integrals b* sm are unity for 
such m values that m, =6;s. 
By comparison with Eq. (74) for the imperfect gas it is seen that the osmotic-pressure equation 
corresponds exactly to that for the pressure of a gaseous solute, except that the activity a, replaces 
t The cluster integrals b*um, in which the solvent subscript is written explicity, should not be confused with the cluster 


- als bnm of the second kind (Eq. (53)) which appear only in Section 7. Also the — subscript » should not be 
used with the chemical potential, ys. 
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the fugacity y,, and the integrals occurring as coefficients are led by the use of the potentials 

of average force in the infinitely dilute solution, involving only the mutual potentials of solute 

molecules and never explicitly the potentials of solvent with solute or solvent with solvent molecules. 
The other equations are: ; 


c.(p*, a) = om m,b*om a,"", (86a) 
p(p*, a) =p* Dim m J], a,”" 
= p* exp B*i,n II. (b) 


1 
In a, =In c,—— n,B*on [], (87) 
Cs 


where the second form of (86b) comes from (81b). 

It is thus to be noted that, for a world at constant temperature at least, the presence of intangible 
molecules which act with forces on the ordinary molecules but which could pass through them, 
would not alter the applicability of the Gibbs statistics to equilibrium systems. Since the intangible 
molecules pass through all barriers their fugacity would everywhere be constant, though not their 
concentration, and all measured pressures would be osmotic pressures. The measured forces between 
ordinary molecules would all be average forces, influenced by the all-pervading fluid molecules, but 
the potentials calculated from them, the potentials of average force, would be strictly applicable to 
the calculation of all observed thermodynamic properties. 

Rather than measuring osmotic pressures, a more usual laboratory practice is to vary the con- 
centrations of solutions under constant (atmospheric) pressure, a process corresponding to the 


equation: 
Ds, m m(Aa)# Il: a,”* = (88) 


The question now arises as to the dependence of one of the variables of Eq. (88) (for example the 
solvent activity increment, Aw) on the other variables (the solute activities), i.e., the determination 
of the coefficients Xi in the equation 


a.". (89) 
(Ai []. 


pi! 


We note that 


(Aa)* (Si wi=n)u! Il 


which, when employed in (88) gives 


Lisui+ms 


i pi! 


im>O i i! 


Lu =0 pi 


where the last two lines are obtained through the substitution k = }°i ivi+m. 
Since (90) must hold for all values of the variable a, the coefficient in brackets must be identically 
zero for each k: 


k 
=0 (all k). (91) 
Mi 


wi 
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(b) 
(c) 
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While for many-component systems this gives fairly complicated relations between the \’s and 
b's for which no general analytic solution has yet been found, the first several \’s are easily eneneneend 
for two components. (91) becomes 


)u: A“*=0 (91’) 


wi 
which for various values.of k yields the relations: 
R=1 AL=—1; 
R=3 og + 11 +20* A2+5* 10A3 =, 


or, 
Ags = (b* 30 — b* 21 12 — b* 03) + (2b* 29 — 11) 11 02 20); etc. 


Thus, for a binary system, the solvent activity as a function of the solute activity at constant pressure 
is given by the expression 


a= p* — b*9)a?+ [(d* 30 — b* 3) + (2b* 29 b* 11) (b* b* o2— b* 29) (92 ) 


The discussion of the application of the general equations to binary systems will be continued in 
the following section. 


10. TWO-COMPONENT SYSTEMS 


Because systems of three or more components have not been so exhaustively studied and are less 
familiar to most than are those of two components, we propose to devote this section to an application 
of the general equations to binary systems. 

In the series notation of Section 8, Eq. (85c), expressing the osmotic pressure in a solution con- 
taining a single solute, assumes the form: 


m(p*,c)=RT[Le— don (n—1)B* onc” | 
=kT[c+G(c, B*on) —cG'(c, B*on) (93) 


where the prime’ indicates differentiation with respect to c. It is seen that for the special case in 
which the solute concentration, c, becomes so-small that powers higher than first be neglected, this 
equation reduces to the familiar law of van’t Hoff. Equation (93) is identical with that obtained 
previously, for example by Mayer and Harrison,* ¢ for the pressure of an imperfect gas, except that 
here the irreducible integrals are evaluated under the conditions which exist in the infinitely dilute 
solution rather than in the infinitely rarefied gas. 

The coefficients, B*,,, involve integrals over the coordinate space of the solute molecules only, 
and the integrands are functionals of the potentials of average force between solute molecules only. 
The functional dependence of the integrands on the potentials of average force of the solute molecules 
in the infinitely dilute solution is exactly the same as the functional dependence of the integrands of 
the virial coefficients of the gas on the potentials between isolated molecules. The presence of the 
solvent does not appear explicitly at all in Eq. (93), and appears only implicitly as the potentials of 
average force of the solute molecules are influenced by the presence of the solvent. 

The equations are derived under the assumption that the integration over the internal coordinates 
of one molecule of the solute converges to a value independent of the volume. This means that the 
assumption is made that dissociation of the solute is negligible. If it were not, an appropriate trans- 
formation of coordinates, changing the description of numbers of molecules, and of “internal coor- 


t Note that here B,,=8,:/n used in this reference. 


f 
é 
t 
t 
\ 


MULTICOMPONENT SYSTEMS 


dinates”’ versus coordinates of center of mass, must be made. In any case, the van’t Hoff law holds 
for the molecular weight appropriate to the solute in the solvent, and holds at low densities as long 
as the potentials of average force between the molecular entities chosen are sufficiently weak to lead 
to small values of the product of ‘“‘virial coefficients’ B*,,, times the solute density (concentration) 
c, measured in molecules per unit volume. Evidently, then, the osmotic pressure for solutions plays 
the role analogous to the total pressure of a gaseous mixture. 

We now examine the conditions under which phase transitions may occur. For our two-component 
system, Eq. (85a) has the form 

w(p*, a)=kT G(a, b* om) =RT b* oma”, (94) 
m21 

where a is the conventional liquid phase activity measured in molecules per unit volume (cf. 
Section 9). 

The slope of the curve obtained by plotting the osmotic pressure against the density of the solute* 
is given by the expression 


an kT ainy 
(=) -— (95) 
Oc/.7 C Olnc/.r 


since a is proportional to y by Eq. (78c). From (86a), 
c= > mb*..0", (96)t 


m21 


Or dln y 
27 


and (95) becomes 


To evaluate the derivative in the right of (97), we note from (96) that 


fi ainy 
OInc/ m>1 OInc/.r 


(=) >, 40". (99) 


m21 


It is now apparent that if the sum in the denominator of (99) has an infinite discontinuity for some 
value g of the fugacity of the solute (corresponding to the concentration 2), the slope of the osmotic 
pressure vs. density curve is zero for all densities above @ which satisfy (96). (The divergence of the 
sum }om>1 mb*,,a" does not necessarily require an infinite discontinuity in the density function 
(96).) If we wish to investigate the slope for fugacities greater than 9, we must use the analytical 
continuation of the divergent sum. 

Returning now to the discussion at the end of Section 8, the singularity, 7, of G(a, b*om) may cor- 
respond either to the vanishing of the determinant in (69), which now becomes 


B*¥ on) =0, (100) 


or a singularity at ¢ in the function G(c, Bon). By (93) we write 


=kT(1—cG’"(c, Bon) 
~0G'"(c, Bon) ] (101) 


* This corresponds to the physical process of adding solute to a given quantity of pure solvent and always readjusting 
the pressure on the solution so as to maintain osmotic equilibrium between the solution and the pure solvent. 
t It is tacitly assumed that the density function (96) increases monotonically in the fugacity, and vice versa. 
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Or 
—) = Box) +06’"(c, Bun) J. (102) 
0c? z,T 


From the shape of the osmotic pressure curve vs. solute density at the phase transitions in a tem- 
perature region far from the critical mixing temperature, particularly its discontinuity in slope, we 
can only conclude that the singularity in the function G(a, bom) corresponds to a singularity in 
G(c, Bon) since otherwise, if € were determined as the root of (100), (@a/dc).,7 would approach zero 
continuously as c approached ¢ from either side. 

In the idealized two-component case for which w,=0 for m>2, the irreducible integrals with 
which we are concerned are defined as 


1 
B Sz II [exp (—taa(z, 0; rer)/RT) —1 (103)" 


(sum over all products with all molecules more than singly connected) 


in which ua(Z, 0; 7,2) is the potential of average force between two molecules of type a (in an in- 
finitely dilute solution of solvent fugacity z) as a function of their distance of separation r;:. The 
presence of the solvent at fugacity z serves only to modify the potential energy function from its 
value at z=0 and not the character of the behavior of the integrals with temperature. 

Although little is known about the behavior of the actual irreducible integrals B,,(z, 0), there 
is no apparent reason why their qualitative behavior with temperature should be radically different 
from that of B,,(0,0), used for the gas. If, as seems probable, this is the case, the arguments of 
Mayer” regarding the nature of the critical region of an imperfect gas may be analogously applied 
to the critical mixing region of two partly miscible liquids. 

Let us then investigate the temperature dependence of certain properties of the osmotic pressure 
function, particularly its first and second derivatives with respect to c, as that variable is increased. 
At the (lower) critical mixing temperature,* T,, defined as the lowest temperature at which the 
first derivative is zero, in order to have a critical point of the ‘‘van der Waals’’ type both first and 
second derivatives (Eqs. (101) and (102)) must be zero. 

If, for all values of T>T., the series cG’’(c, B.,) approaches a value less than unity as c approaches 
@, then ¢(7) will be the density at which a second phase appears, and T,, will be the (lower) limiting 
temperature at which a second phase may be obtained. If it happens that 


lim Bon(T-) ]<1, (104) 


the slope given by (101) will not be zero at the critical mixing point—a type of critical point 
which is never observed. Furthermore, since it seems unlikely that the values of the B,,’s for a 
particular system should be such that at the critical mixing point we should have fulfilled simulta- 
neously both the conditions (for c=é) : 


cG"'(c, Bon) =1 (105) 


and 
CG" (c, Bon) = —1, (106) 


which presumably may be accomplished only by very special types of attractive forces, it appears 
highly improbable that al/ critical mixing points are of this type. For the case in which the limit (104) 


17 Reference 1, pp. 308-314; reference 10, pp. 639-640. 

*In liquid systems some confusion may arise due to the fact that both ‘‘upper’’ and “lower” critical mixing points 
may occur, and even in the same system. In what follows we shall consistently speak of the lower critical point, but what 
is said will be equally true for the upper critical point with only slight changes, such as the order of certain temperatures: 
T-<Tm (lower); T->Tm (upper). The subscripts c and m are those used in reference 1, and should not be confused with 
concentration c, which never occurs as a subscript, or the running index m. . 
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is greater than unity the Jacobian (69), changes sign; the expansion (94) is no longer valid, and Eq. 
(101) would indicate a (physically ridiculous) osmotic pressure which decreases with increasing solute 
density. The only other apparent alternative is that the sum in (105) becomes unity at some value 
of the temperature, 7;,, which is higher than the true critical temperature, and that for temperatures 
between 7. and T,, the singularity in G(a, bo») is determined by condition (100) rather than by a 
singularity in G(c, Bo»). The solution of (105) satisfies the condition (106) at T=7.. An anomalous 
first-order transition,'® in which the r—c curve is horizontal over a finite range of density but has 
no discontinuity in slope, thus occurs between the temperatures 7,, and T,. An analysis strictly 
analogous to that used by Mayer for the imperfect gas'® shows that, under the hypothesis that the 
B.» have the same sort of temperature dependence as the 6;, the separation into two conjugate 
phases occurs only at temperatures greater than T;,,, at which point the surface tensions of the two 
phases are identical. The argument will not be repeated here, but it should be noted that it is the 
mutual surface tension between the two phases which should be considered. 

Probably because of the extremely high osmotic pressures involved in all but the most dilute 
solutions, there are no such measurements of osmotic pressure with which to compare the foregoing 
theoretical description of the region of the critical mixing point. However, constancy of the osmotic 
pressure requires constancy of the partial vapor pressures of the two components over the solution, 
since:the vapor pressure, p, of the solvent in the solution is uniquely expressed in terms of the 
osmotic pressure as the solution of the equation: 


Po B B ‘ 


Pp 


in which v, is the molar volume, p, the vapor pressure, and 8 the compressibility of the pure solvent; 
and 6(p) is the partial molar volume of the solvent in the vapor phase above the solution. The region 
of anomalous first-order transition is thus characterized by a finite concentration range in which 
the plot of the partial vapor pressure of either component (or, for that matter, the total pressure of 
the solution) against concentration is horizontal but has no discontinuities in slope. 

There is much experimental evidence supporting this conclusion.* As early as 1900 van der Lee*! 
remarked on the flatness of these curves for the system: phenol-water. The same characteristic of 
the system: cyclohexane-aniline was noted and briefly discussed by Mme. H. Schlegel.” Measure- 
ments by Dr. L. D. Roberts” in this laboratory showed that over a wide range of composition the 

‘partial pressure of triethylamine in triethylamine-water mixtures, at temperatures a few degrees 
below the (lower) critical mixing temperature, was constant to within approximately 0.1 percent. 


11. SUMMARY 


The method of the grand-canonical ensemble has been generalized for multicomponent systems, 
and equations have been obtained for the calculation of the thermodynamic properties and dis- 
tribution functions of such systems at any particular set of fugacities in terms of those functions 
at any other fugacity set. The theory is thus applicable to liquid solutions, for which the osmotic 
pressure has been shown to play a role analogous to that of the total pressure of a gaseous system. 
It has further been shown that the existence of an anomalous first-order transition near the critical 
mixing temperature of a binary mixture is probable. 


18 Cf. J. E. Mayer and S. F. Streeter, J. Chem. Phys. 7, 1019 (1939). 
19 See reference 1, pp. 310 ff; reference 5, pp. 93-94. 
20 Cf, for example, A. A. Noyes and M. S. Sherrill, Chemical Principles (The Macmillan Company, New York, 1938), 
second edition, paragraph IV, 25. 

* The anomalous behavior of imperfect gases, with experimental comparisons, in the critical region has been thoroughly 
discussed by S. F. Harrison and J. E. Mayer, J. Chem. Phys. 6, 101 (1938). 
*t Zeits. f. physik. Chemie 33, 628 (1900). 
2 H. Schlegel, J. Chem. Phys. 32, 215 (1935). 
*3L. D. Roberts and J. E. Mayer, J. Chem. Phys. 9, 852 (1941). 
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APPENDICES 
A. Thermodynamic Excursion 


Since the independent variable z,, which we call fugacity,* is not a common one, it may be well 
to discuss briefly its properties. As in Eq. (18), we may write 

Ms =kT In (A1) 

where yos(7), a function of the temperature alone, is the chemical potential of a perfect gas of pure 


constituent s at a density of one molecule per unit volume. The total derivative of u..(7) with 
respect to T may be calculated as ' 


= p+ pure s, =1 (A2) 


= —Sotk, 
where’ S,, is the entropy per molecule of pure perfect gas of species s at unit density, and k is the 
Boltzmann constant, the latter following from the fact that the perfect gas equations are to be em- 


ployed for pos. 


From (A1) we then have 
(A3) 


and 
In .+dyo,/dT 


1 
= Mos t+kT — 


1 
T'S ]p,=1, pure gas} 


1 


where v is the volume per molecule and £,, is the energy of a perfect gas of molecules of type s at 
unit density, which is, of course, a function of temperature alone. We also have that yu, is a function 


of z, and T alone, or, 
(Ous/02,)T,2=0, (AS) 


Let us now consider the extensive thermodynamic potential : 


PV=F-A=). Nu.—A, (A6) 


d(PV)=D. Nedust wd N.— udN,+PdV+SdT 
=>. N.du.+PdV+SdT, 
from which it is seen that PV is a “natural” function of u, V, and 7, namely that 
[O(PV)/du. |v, Tur=N,, (A7) 
[a(PV)/dV]r,.=P (A8) 
[a(PV)/dT lv, (A9) 


Corresponding to this, the grand partition function for a grand ensemble of fixed V, T, and yw is 
equal to the exponential of PV/kT (Section 4). 


* The fugacity of Lewis and Randall is proportional to z,, the former becoming equal to the pressure at infinite dilution 
instead of to the number of molecules per unit volume. The Lewis and Randall fugacity is thus kT2s. 
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From (A7, 8, and 9) it is seen that the intensive property P =(PV)/V, if expressed as a function 
of the intensive variables 7, w has the partial derivatives 
(0@P/dT),=S/V, (A10) 
(OP /Ous)T, ur = N;s/V=ps, (A11) 


and, of course, (@P/dV)7r,,=0. 
Now let us examine P as a function of 7 and z, using (A4, 5, and 10) to obtain 


Les (OP/ Ape) 


1 
Ls Ne(us— Eos) ] 


=(TS+F-E,)/TV 
=(H-E,)/TV =(PV+E-E,)/TV, (A12) 


where E,=>_.N.E,s is the energy of the constituent molecules of the system at unit density, and 
H=E+PYV is the enthalpy. Using (A3, 5, and 11) one arrives at 


We may regard 2, as a function of the intensive variables P, T, and z, for all r#s. By inversion of 
(A13) we find 


(0 In 2,/0P)7, z-=1/kT p, (A14) 


while the negative ratio of (A12) to (A13) gives 
(0 In 2,/0T)P, 2-= 


(A15) 


and the negative ratio of two equations of the type (A13) shows that 


(0 In In 2,)7,P,2-= — pr/ — (A16) 


The second derivative (d°P/dT0Z,)z, may be calculated by differentiating (A12) with respect to 
zs, or (A13) with respect to 7. Equating the two results leads to the equation 


(Op./8T )z= (2./ VRT*) [AE (z) /02, ]7, zr. (A17) 


If 2, is regarded as a function of P, T, and N,, Eq. (A1) permits us to write all derivatives from 
common forms, leading to the familiar Gibbs-Duhem relationship for all changes at constant T and P: 


N.d In pd In 2:)r,p=0, (A18) 
showing that z, depends only on N,/N, of which Eq. (A16) is seen to be a special case; furthermore, 
(0 In 2,/0P)7, N,=d,/kT, (A19) 


where 3,=(0V/0N,)P, 7,N. is the partial molecular volume. The derivative with respect to tem- 
perature requires a little computation, and is found to be 


(0 In Mos)/RT? 
= —(A,—Ho.)/kT? (A20) 


where H/, is the partial molecular enthalpy of species s, 1,=(0H/AN,)r,7,n,, and H,, is the en- 
thalpy per molecule of a perfect gas composed of molecules of species s at unity density. 
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B. Reversion of the Expansion of the Potentials of Average Force 
We wish to solve the set of equations 
Wn(z, {N})=L wn(z, {n}N), N=1, 2, (B1) 


discussed in Section 6 (Eq. (47)) for the component potentials of average force wn(z, {n}). Let us 
first note that each molecule is assigned a number, so that they are all distinguishable ; if we associate 
with each molecule a second number denoting species it makes them none the more distinguishable 
(except, of course, when we perform integrations, in which case coordinate interchanges are limited 
to pairs of molecules of the same species). Therefore, for the purposes of this and the following section 
we may simply consider a set of numbered molecules without regard to species, and not write the 
symbols in bold-face; also, since all the quantities are functions of z, this will be omitted for con- 
venience in writing. 
Let us consider, then, the set of equations 


= D2 Jnwet +L (B2) 


W,(1)=w,(1), (BS) 


where the symbol >-[—z7], indicates summation over only those subsets ®f » which are n—7 in 
size. Evidently there will be one equation of type (B2); n-equations of type (B3), depending upon 
the choice of the unique molecule included in the set {v} but not in the subset {n—1},; and in 
general, n!/(n—1) ¢!-equations of type (B4). 

If we were to sum over all the n-Eqs. (B3) and subtract (B2), the terms w,-; would cancel. We 
propose now to generalize this procedure so that all terms on the right cancel except wy. 

In the set of Eqs. (B4), the number of equations having the term w,_; for {n—j},-; particular 
molecules (of the set »—7) is just the number of ways in which we can choose i-molecules from the 
set {m} without including any of the set {n—j}—or, the number of ways in which we can choose 
the set {7} of molecules from the set {7}; this is simply j!/(j—7)!7!. Again employing the symbol 
>[n—j]n to indicate summation over only those subsets of which are n—j in size, 


7=i 
Now we form the sum 


n—1 


stated n-12>j2i21 


n—1 j 


j=1 i=1 


n—1 


=> (1-1)! 


wa(z, Wwe, {N}a), 


which is Eq. 
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B7) 


B8) 
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C. Reversion of the Expansion of the Distribution Functions 
In Section 7 (Eq. (49)) we wrote 

k 
F(z, (M})=2 TT gm(z, {m;}M). (C1) 


As in Appendix B we may omit the variable z and the bold-face without loss of generality. This 
device aids the understanding of the proof as well as the explanation. 

By writing out the series (C1) for M =1, to 5, it is possible to guess the solution for the functions 
g, in the form 


l 
=D 0-12) TT Fut M},. (C2) 


j=1 


This step is straightforward but lengthy, and we shall not include it here. We shall prove the cor- 
rectness of (C2) by substitution with (C1) and reduction to an identity. This substitution gives 


U k 
= u(—)U—1)! IT Ms} II toms} 


in which P/* is the number of ways in which K-particular subsets {m;} may be grouped together in 
exactly /-sets—or, the number of ways of putting K-distinguishable objects into /-numbered boxes 
with at least one to each box. We have found no closed form for the general quantity P;*, but the 
following recursion formula, valid for K>2, is easily obtained : 


1 


4... 
PA=P,,+I/P: , (C4) 


in which the first term on the right corresponds to putting the Kth-subset alone in the /th-box, and 
the second term, to putting it into any one of the / boxes which already contain K—1 subsets. 
Now the coefficient in (C3) becomes (for K 2 2) 


K K K-1 K-1 


l=1 l=1 


K-1 K-1 
=> 


=0 (K22), 


K-1 
= > 
l=1 


since P;* is zero for /=0 or 1>K. P;' is obviously unity, and consequently the coefficient is given by 


(—)10-1) = (C5) 


l=1 


and Eq. (C3) thus reduces to an identity, since for K=1, {m;} =. Reverting to the bold-face 
notation, we obtain 


gm(z, 4(—)1(—1)! J] {M,}m), (C6) 


i=1 


which is Eq. (50). 
The inverse, (C2), has been previously reported’ but this appears to be the first published proof. 
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D. Glossary 


For easy reference the following compilation 
of symbols is appended (the numbers in paren- 
~ theses refer to the defining equations) : 


a 


Sais 


F 
Fa(z, {n}) 


gm 


gnm 

G, 

h 

1 

1, j, k,l 
In 


J{x/q] 


k 
k, K,l 


Activity of the solute in a two-com- 
ponent system (94). 


Activity of the solute species s (78c). | 


Helmholtz free energy (Appendix A). 
The cluster integrals of the first kind 
(52); bo =0. 
Cluster integral evaluated at zero 
fugacities (76). 

Cluster integral evaluated at fugacity 
set z* (80, ff.). 

Cluster integral of the second kind 
(53); employed in Section 7 only. 
Concentration of the solute in a two- 
component system (Section 10). 
Concentration of solute species s 
(Section 9). 

Thermodynamic ‘internal energy”’ 
(Appendix A). 


The energy of the /th quantum state i 


of a canonical system. 

The number of degrees of freedom 
(39). 

The force on molecule 7, along coor- 
dinate g (Section 6). 

The Gibbs free energy (Appendix A). 
The distribution function (Section 
2b). 

Correction term for m_ proximate 
molecules (49, 50). 

Second type of correction term (51). 
Defined by (62). 

Planck’s constant. 

“Internal,”’ as on Q;, (1”). 

Running indices. 

An auxiliary function (Section 5), 
(36). 

The Jacobian between the coordinate 
frames g and x (11). 

Boltzmann’s constant. 

The total number of subsets for a 
particular subdivision (49, 50, C3) 


m, M,n, N Sets of molecules. indicates 


summation over all sets for m2 0. 
The spatial coordinates of M mole- 
cules. 


m, M 

P, P(z) 
PONV 
PF 

{Pn} 


Xisy Visy is 


Numbers of systems (Section 4). 
Pressure. 

Defined in Eqs. (25), (28). 

Appendix C. 

The momenta of N molecules (39). 
Coordinate g of molecule 7, (42), (7). 
Canonical partition function (Eq. (3) 
and Section 4). 

Single particle partition function (1). 
Single particle internal partition func- 
tion (1’’). 

Distance between the centers of mass 
of two particles (Section 6). 
Designation of species of molecules. 
(Subscripts.) 

Entropy (Section 4, Appendix A). 
Absolute temperature. 

Component potential energy (Sec- 
tion 6). 

Designates “unconnected” subsets 
(Eq. 49 etc.). 

The potential energy (Section 6), 
(40). 

Volume. 

Partial molar or molecular volume 
(A18), (107). 

Component potentials of average 
force (Eqs. (47), (48)). 

The potential of average force (Sec- 
tion 6). 

Cartesian coordinates of the center 
of mass (Section 2a), (7). 

The fugacity or fugacity increment 
of component s. 

Equation (67). 

The fugacity of component s. 
Equation (66). 


Solvent activity (Section 9). 

The irreducible integrals used in ref- 
erence 5. 

The irreducible integrals, implicitly 
defined by (60); Bo=0. 

The irreducible integral (Section 9) ; 
Boo =0. 

The activity coefficient (19). 

Solvent activity coefficient (77d). 
Kronecker’s delta. =0 for i#j7 and 
=1 for i=). 
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The energy of a single molecule (Sec- 
tion 1). 

Set or subset of molecules, Appendix 
C; also a running index, Section 9. 
A chemical potential set (16); it 
never appears in brackets as inde- 
pendent variable as does the subset 
u of molecules. 

The chemical potential of species s 
at density po. (18). 

Indicates a product over all s. 

A density set (14). 

Solvent density (Sections 9, 10). 
The total number of components. 
Summation. 

A volume element in configuration 
space (Eq. 103). 

The eigenfunction corresponding to 
quantum state / (Section 4). 


w, Q 
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Number of quantum states (Section 
4). 


Special symbols 


A bold-face symbol indicates a set, one for each 
species (Eq. (9)). 


{n} 


new 
{k{mjm}. 


The coordinates of the set n of 
molecules (Eq. (10)). 

A kind of dot product (Eq. (20)). 

A set of k unconnected subsets {m} 
of the set M of molecules. 

oth-order determinant (Section 8). 
Average value of (Section 6). 
Superscript, refers to the standard 
state, z* (Section 9). 

Superscript, refers to the state z=0 
(Section 9). 
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Letters to the Editor 


HIS section will accept reports of new work, provided 
these are terse and contain few figures, and especially 
few halftone cuts. The Editorial Board will not hold itself 
responsible for opinions expressed by the correspondents. 
Contributions to this section should not exceed 600 words in 
length and must reach the office of the Managing Editor not 
later than the 15th of the month preceding that of the issue in 
which the letter is to appear. No proof will be sent to the 
authors. The usual publication charge ($3.00 per page) will 
not be made and no reprints will be furnished free. 


The Fluorescence of Naphthacene in Anthracene 


E. J. BowEN 
Physical Chemistry Laboratory, Oxford, England 
June 4, 1945 


N a recent communication to this Journal,! Dr. Ganguly 
has questioned the need to assume an ‘‘exciton” process 
in the fluorescence of solid solutions of anthracene and 
naphthacene. His experiments show that the dissolved 
naphthacene fluoresces when excited by blue light which 
it absorbs. This is clearly a case of direct excitation of the 
naphthacene molecules. When an exciting wave-length of 
3650A is used, however, the phenomena are not so simple. 
Pure anthracene fluoresces in this light a violet color 
(bands 4030—4450A). A solid solution with as little as 
1 mole naphthacene to 10‘ moles anthracene however 
fluoresces with the characteristic green color due to the 
naphthacene, and at a mole ratio of 5X10‘ the green is 
very strong and the anthracene emission largely suppressed. 
There are only three possibilities for this behavior: 

(1) The naphthacene absorbs more light than the 
anthracene at 3650A. This is ruled out by the fact that the 
extinction coefficient of anthracene is several times larger 
than that of naphthacene at this wave-length, so that the 
ratio of absorptions is about 10*:1 in favor of anthracene 
for the above solid solutions. 

(2) The anthracene absorbs the light and emits its 
violet fluorescence which is taken up by the naphthacene 
and re-emitted as green fluorescence. This is likely to 
occur only when the absorption of the anthracene fluores- 
cence by the naphthacene is at least as great as the absorp- 
tion of 3650A by the anthracene crystal. Using extinction 
coefficient data obtained for liquid solutions of these 
hydrocarbons, a maximum value for the anthracene-naph- 
thacene ratio satisfying the above condition is about 10, 
or one thousandth of the actual ratio. 

(3) The anthracene absorbs the light and hands the 
quanta of excitation energy on from molecule to molecule 
through the crystal by a process of resonance until either 
an anthracene molecule loses it again by fluorescence or 
until it is trapped by a naphthacene molecule which is 
thereby caused to fluoresce. 

The latter process explains why crystalline anthracene 
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has a half-life for fluorescence of about 10-5 sec. while a 
liquid solution has a half-life of 10-® sec. These figures 
indicate that the distance the excitation energy (exciton) 
travels between absorption and emission for pure anthra- 
cene is of the order of 1000 molecules. At naphthacene 
contents of 10~* there would be one molecule of naph- 
thacene for every 30 of anthracene in any linear direction 
in the crystal so that even at this low concentration there 
is a high chance of the energy reaching a naphthacene 
molecule. 

Solution of these crystals in a liquid destroys the 
“exciton” effect, which is of course dependent on proper 
orientation of the molecules. A strongly green-fluorescing 
solid solution, when dissolved in benzene, shows only the 
blue fluorescence of pure anthracene. The fluorescence of 
solid organic compounds are liable to be much more sensi- 
tive to traces of impurities.embedded in the crystal than 
the corresponding solutions in liquids. The hydrocarbon 
coronene, for example, fluoresces blue in solution, but even 
highly purified crystals extracted from tar fluoresce green ; 
whence it is probable that these crystals still contain a 
minute quantity of impurity not detectable by ordinary 
organic chemical means. 

The brilliance of the green fluorescence of anthracene 
containing a few parts in a thousand of naphthacene 
renders the substance very suitable for preparing fluores- 
cent screens for focusing U.V., spectrographs. A hot 
acetone solution of the hydrocarbons is added to about 
three times its volume of warm dilute gelatine solution in 
water and the mixture allowed to dry in not too thin a 
layer on a glass plate. Such plates are translucent and 
permit sharp focusing. They are also useful for rapid visual 
observations with a U.V. spectrograph. 


1S. C. Ganguly, J. Chem. Phys. 13, 128 (1945). 


The Ratio Rule 


WALTER F. EDGELL 
Department of Chemistry, State University of Iowa, Iowa City, Iowa 
May 17, 1945 


OETHER! has recently found that the results of his 
studies on the isotopic methyl halides could be 
expressed very well by the following empirical ‘ratio 


rule:”’ cp.x 
Vn 
(1) 


where X is a halogen atom and », refers to the corre- 
sponding frequency in each molecule. Moreover, he has 
shown that this rule could be successfully extended to 
other pairs of isotopic molecules as long as they were of 
similar symmetry. No theoretical derivation of this rule 
was known, although the results suggest that here is a 
phenomenon of considerable importance to the estimation 
of the effect of isotope substitution on molecular vibrations. 
It is the purpose of this note to indicate briefly the 
theoretical background of the ratio rule and to propose, 


ynCD3Cl = 
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in addition, a more useful but generally less accurate 
form of it. 

The matrix expression for the frequencies of vibration 
of a molecular system is given by? 


\=LGFL, (2) 


with \ a diagonal matrix whose elements are \ii=42v;?; 
L is the transformation matrix from symmetry coordinates 
to normal coordinates, and G and F are respectively the 
kinetic and potential energy matrices. This equation may 
be regarded as a set of relationships between each \;; and 
the various L;;. If the equation involving dj; is differ- 
entiated partially with respect to L;:, one obtains after 
some manipulation 


This is Rayleigh’s principle and shows that \,; is stationary 
for small variations in the normal coordinates. Thus 
when a somewhat approximate normal coordinate is 
available, Eq. (2), or its equivalent, yields a rather accurate 
value for the corresponding frequency of vibration. 

These results may be applied directly to the isotope 
effect by taking each normal coordinate of the non-isotopic 
molecule as a first approximation to the corresponding 
normal coordinate of the isotopic molecule. Using Lo, \° 
etc., and A, G etc., to refer to non-isotopic and isotopic 
matrices respectively, one obtains from (2) 


N= (Lo GL 
= (4) 
Equating corresponding diagonal elements yields 
= KimKinGmn®, (5) 


xix being most simply thought of as the unnormalized coeffi- 
cient (Lo). A more formal perturbation treatment 
(similar to that once proposed by Wilson‘) shows that 
this equation is correct to the first order. 

When symmetry coordinates are formed from internal 
coordinates (bond stretches, angle deformations, etc.), 
they can be made to correspond rather closely to normal 
coordinates. If the symmetry coordinates of the non- 
isotopic molecule have been so chosen, the principal terms 
in Eq. (5) involve «i:2; and one may write that 

KiPGii/ 
=Gi/Gii®. (6) 

Equation (5) is the fundamentally basic equation of the 
ratio rule (correct to the first order). Equation (6) is un- 
doubtedly the most useful form of the rule, its field of 
application being limited only by the assumptions used in 
deriving it. Thus it is equally valid in calculating the 
frequency shifts resulting from the change CCl,* to 
CC1;CI8? as for CH, to CD, both of which fall outside the 
scope of Eq. (1). 

On the other hand Eq. (1) is the form of Eq. (5) which 
is valid whenever two pairs of molecules have essentially 
the same normal coordinates. Unfortunately this require- 
ment restricts the use of Eq. (1). However, when these 
conditions are met, as with the two pairs CH;Cl, CD;Cl 
and CH;Br, CD;Br, Eq. (1) yields more accurate results 
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than Eq. (6). This, of course, is the consequence of the 
fact that symmetry coordinates are not quite normal 
coordinates, a fact more fully corrected for in Eq. (1) 
than (6). 

The kind of accuracy that may be expected from Eq. (6) 
under favorable conditions is indicated by Table VIII of the 
Noether paper, in which he notes the agreement between 
the frequency ratio and the square root of the G ratio for 
the methyl halides. However, calculations show that Eq. 
(6) does not give good results when a D-H substitution is 
accompanied by a symmetry change (i.e., CH -CH;D). 
This is the result of near degeneracy in the A; and E sym- 
metry classes; an approximate perturbation treatment 
appropriate thereto gives excellent results. It should also 
be pointed out that whenever a symmetry change is in- 
volved special symmetry coordinates should be used in 
setting up the G matrix. This whole topic will be discussed 
in greater detail in a forthcoming paper on this subject. 

Since the ratio rule deals with individual frequencies of 
vibration, it may be expected that it will be even more 
useful than the well known Teller-Redlich product rule.® 


1H. F. Noether, J. Chem. Phys. 11, 97 (1943). 

2 E. B. Wilson, Jr., J. Chem. Phys. 9, 76 oc 

3See R. A. Frazier, W. J. Duncan, A. R. Collar, Elementary Matrices 
(Cambridge’ University Press, New York, 1938), p. 299, for a closely 
related discussion involving ‘‘modal columns.” 

. Wilson, i Phys. Rev. 45, 427 ag 

5 Teller, e. <. K. Ingold et al., J. Chem , London, 971 (1936); 
Zeits. f. physik. Chemie B28, 371 E. B. Wilson, Jr., 
reference 


Collisional Deactivation of Excited 
Ethylene Molecules 


DAPHNE LINE AND D. J. LERoy 
Department of Chemistry, University of Toronto, Toronto, Ontario 
May 28, 1945 


N the reaction of Hg(*P:) atoms with ethylene it has 

been shown"? that acetylene and hydrogen formation 
can be explained on the basis of the reactions 

+Hg('So), (1) 

(2) 

(3) 

where C,H,* represents an excited ethylene molecule. This 

mechanism was based on the following experimental facts: 


(a) oa and hydrogen formation is not inhibited by nitric 


oxi 
(b) suling allowance for the change in collision frequency with 
temperature for the bimolecular reaction (2), the rate of 
formation is independent of temperature, and 
(c) the rate decreases with increasing ethylene pressure. 


At 25°C almost all of the ethylene used up in the early 
stages of the reaction is accounted for by the acetylene 
and hydrogen formed. This makes it relatively easy to 
investigate the deactivation process, (2), at that tem- 
perature. If this reaction occurs then qualitatively at least, 
it is immaterial whether the deactivator is an ethylene 
molecule or not. Carbon dioxide was chosen as a deactivator 
because of its small quenching cross section relative to 
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ethylene (3.5A? vs. 48A*) and the fact that quenching by 
carbon dioxide does not lead to decomposition.’ 

A series of runs was first made to locate the maximum 
rate.! This was found to occur at 4-mm ethylene pressure, 
the average rate at this pressure being 0.012; mm per min. 
A number of runs with 4 mm of ethylene and 12 mm of 
carbon dioxide gave an average rate of 0.006;. From the 
quenching cross sections of carbon dioxide and ethylene it 
follows that under these conditions only 85 percent of the 
Hg(?P:) atoms collide with ethylene molecules and since 
it was found that the rate was directly proportional to 
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the light intensity the rate in the presence of carbon 
dioxide must be compared with 85 percent of 0.0127 or 
0.010;. We thus have further evidence to substantiate 
reactions (1)—(3). Furthermore, a comparison of the rate 
in the presence of carbon dioxide with 85 percent of the 
rate with 16 mm of ethylene alone (0.006,;) shows that 
carbon dioxide is as efficient as ethylene in deactivating 
excited ethylene molecules. 


1D. J. LeRoy and E. W. R. Steacie, J. Chem. Phys. 9, 829 (1941). 
2D. J. LeRoy and E. W. R. Steacie, J. Chem. Phys. 10, 676 (1942). 
3 J. E. Cline and G. S. Forbes, J. Am. Chem. Soc. 61, 716 (1939). 


